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AUTOMATIC CONTROL AND THE IRE 


This is the first Transactions of the Professional Group on Automatic Control. 
Like the Group it is small, and its infancy is reflected within by the type- 
written, offset printing. When the Group grows larger, more mature and in- 
creasingly active, the Transactions will be typeset and more attractive. 


Although the PGAC is young, the IRE actually has been active in the field 
of automatic control for a number of years. In the past some of the first 
_ theories of feedback control were developed by those interested in relatively 
simple electronic feedback amplifiers. Radio and electronics engineers were 
among those who created new theories and developed more advanced feedback cir- 
cuits until today radios, radar, television and computers contain numerous 
complex feedback control circuits which include semiconductors, electromagnetic 
devices and mechanical elements as well as vacuum tubes. 


And a new trend is developing -- all of the complex marvels, the radio, 
the radar, the television, the computer, each containing internal feedback 
control loops, are being combined into over-all automatic control systems. 
These have been referred to by the military in the weapon systems concept and 
by industry as automation. Unique problems have appeared during the integra- 
tion of these complex systems; new theories such as those concerning the 
nature of information which must flow between system elements, and new tech- 
niques such as those involved in digital computation are being developed. 


In the future automatic control will become more complex; it will encom- 
pass broader fields. There will be a greater need for developing and incor- 
porating new theories, techniques and components into integrated systems. The 
IRE may well participate in the evolution of these integrated systems, not only 
because of the past association with feedback control but because at the pres- 
ent time many of the new techniques and theories necessary for the future de- 
velopment of automatic control are being fostered by other professional groups 
in the IRE. It will be a function of the future automatic control engineer to 
pose new problems for specialists in the related groups and to coordinate new 
information into integrated, automatically controlled systems. 


To represent the automatic control engineer, to provide a means for dis- 
seminating information about automatic control, its related interests, its 
problems and its developments, the Professional Group on Automatic Control was 
formed. From the constitution of the PGAC: "...The field of interest of the 
Group shall be automatic control systems. It shall encompass the components 
thereof, such as transducers, data transmission links, computers and control 
devices and the integration of these components into control systems...." 


Although the interest of the Group will concern the problems of combining 
elements into over-all complex systems and the solutions of these problems, at- 
tention will be given to any subject, abstract or practical, relating to auto- 
matic control within any of the system elements, simple or complex. It should 
be observed that.the field of interest is broad and the subjects to be consid- 
ered varied. This first Transactions issue illustrates the range of interest. 
Although there are only seven papers, the topics covered are diversified, from 
the first which is a specific example of the integration of a digital computer 
into a feedback control system, to the last which is a general discussion of 


the future of automatic control and its relationship to the engineer. 


The intervening papers concern linear and nonlinear theory applied to 
Simple feedback systems and to practical solutions of problems inherent in 
the development of a complex fire control system, 


Each of the papers was selected from a large number by a Papers. Study and 
Procurement Committee. Although most of these papers were conference papers, 
not specifically intended for the Transactions, others are being obtained from 
local PGAC chapters where they may or may not have been formally presented. 
Actually, any papers concerning automatic control will be welcome, but they 
will have to be reviewed by the committee before being published in the Trans- 
actions. Rejected manuscripts will be returned to the author with comments 
and suggested revisions. To aid in preparing papers and to minimize rejects, 
a Suggested form for writing manuscripts is included in this issue. Papers 
should be submitted to Mr. M. R. Aaron, Bell Telephone Laboratories, Murray 
Hill, New Jersey. 


The growth of the PGAC and the aid which the IRE can give to the develop- 
ment of automatic control will depend considerably on the activity of local 
PGAC chapters and the dissemination of information within the Group and among 
other groups. The Transactions is the medium whereby this information can be 
distributed; it is the agent where new problems and developments can be pro=- 
posed. Its usefulness will depend on the material it will contain: articles 
on general theory and new practical applications, tutorial articles, historical 
surveys, design methods, discussions of papers already presented, abstracts of 
papers printed in other group Transactions which may be of interest to PGAC 
members, book reviews, chapter activities or theoretical problems for educa- 
tional purposes. Comments and suggestions for improving the Transactions in 
any way, by changing the format or in the selection of material, will be most 
welcome == and necessary -- if the Transactions is to serve its purpose in the 
development of automatic control. 


SYNTHESIS OF A NONLINEAR CONTROL SYSTEM*¢ 


I. Flugge-Lotz and C. F. Taylor 
Stanford University 
Stanford, California 


Summ -- The investigated nonlinear control system consists of a linear mem- 
ber with an ensemble of possible discrete combinations of proportional and 
derivative feedback around the linear member. The particular combination of 
proportional and derivative feedback employed at any instant is determined by 
a feedback switching circuit which is in turn operated by sensed binary infor- 
mation obtained from the output, output derivative, error and error derivative 

namely, the signs (sgn) of these variables]. Techniques that are common to 
the digital computer field are used to implement this switching circuit. 


For a linear member of second order the feedback circuit is comprised of 
four discrete values of proportional feedback and four discrete values of 
derivative feedback. 


Simulation techniques have been used to study and evaluate the perform- 
‘ance of the nonlinear control system and to compare it with a linear system 
for a wide variety of inputs. A sample of these experimental results is pre- 
sented. The experiments show that the nonlinear system performance is much 
better than that of a linear system of comparable power handling capability. 


Preliminary studies of a nonlinear control system of third order show 
that the basic idea can be extended to higher order systems. 


The decision to introduce nonlinear components in a control system may 
stem from the desire to have a system whose performance is less tied to his- 
tory than that of linear control systems of comparable power handling 
capability. 


Such a nonlinear system has been mathematically described earlier.¢ In 
the present paper its physical realization and experiments for studying per- 
formance are described. 


The particular. control process may be visualized as 


(1) establishing an ensemble of linear differential equations with constant 
coefficients (some with stable and some with unstable homogeneous solu- 
tions) and 


(2) switching from one member of the ensemble to another so as to maintain 
small instantaneous error between input and output for relatively 
arbitrary inputs. The switching sequence and time of switching are 
functionally dependent upon quantized binary information of the output 


Research supported by National Advisory Committee for Aeronautics (NACA). 

is paper was presented at WESCON, the West Coast Conference sponsored by 
the IRE - West Coast Electronics Manufacturer's Association, in August, 1955. 
$ritigge-Lotz, I., and Wunch, W. S.: "On a nonlinear transfer system." Jour. 


Appl. Phys., pp. 48-488; Apr., 1955. 


and input, namely the signs (sgn) of the output, output derivative, 
error and error derivative.’ 


Briefly, the physical manifestation of this control process is a linear system 
with discrete variable proportional and derivative feedback around the linear 
member. The particular combination of proportional and derivative feedback 
employed at any instant is determined by a binary logic switching circuit. 
This circuit operates with sensed quantized information of the output and 
input. 


MATHEMATICAL DESCRIPTION OF A SECOND ORDER NONLINEAR CONTROL SYSTEM 


Mathematically, the system to be considered is described by 


D2 
— + 2D(1 +6) =. + (1 +%)y = x(t) (1) 


where x(t) is the input, y(t) is the output and D is the linear damping factor 
(when ba aL AE ie 0). The instantaneous error is defined ase =x -y. The 
parameters B,, and Y, are stepwise switching functions of their implicit vari- 
able, time, and are given explicitly as functions of x and y. That is, 


Bn = -sgn(y*) ig sgn(e) ne 26 sgn(e')] 


63 = 38+ of 
8» a 6 = ob ,6 and of are 
6, = -Py 26 = as positive constants 


Bo = 8 - 26 = -6, 


Yn = ~sen(y) [yy sen(e) + zy sen(e*)] 
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PHYSICAL REALIZATION 
In order to understand the nature of a physical: system described by Eq. 


(1), this equation is rearranged and rewritten in operational form (utilizing 
the Heaviside operator p = d/dt) as 


[p? + 2p +1] y 2x - [e0Ap +ynly | (2) 
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From Eq. (2) a symbolic representation of the nonlinear system is obtained in 
the block diagram form as shown in Fig. l. 


With the aid of the block diagram the interpretation of Eq. (1) is 
straightforward. The nonlinear system consists of 


(1) a linear second order member. 


This could be a simple position servo, 
for example. 


(2) a feedback circuit comprised of 


(a) four discrete values of proportional feedback (two positive and two 
negative) denoted by gain constants,y,. 


(b) four discrete values of derivative feedback (two positive and two 
negative) denoted by gain constants,4,. 


The particular (Bns Vn) feedback combination employed at any instant is 
determined by a switching logic circuit. 


SWITCHING LOGIC 


In the physical realization of Eq. (1) a switching logic has to be estab- 
lished that ensures that the proper (Oe, Yn) combination is employed at the 
proper time. By definition, the parameters 8, andy, are functionally depend- 
ent upon the signs (sgn) of the four variables y, y', e and e'. Since the 
sign of a variable is binary information it is convenient to utilize digital 


p* tacdp #1 OUTPUT 


SWITCHING 
LOGIC 


Fig. 1 - Block diagram of Eq. (1). 
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computer techniques in establishing the necessary switching logic. This is 
done as follows: 


Let y > 0 be represented by "0" (binary zero), and 
y ¢ 0 be represented by "1" (binary one). 
Similarly for y*, e and e*. 
If the ordered sequence is established as 
[y, a Aled e4] 


then the allowed (6,,, y,) combinations may be encoded in four digit binary 
logic (namely, binary coded decimal). This is done in Table I. 


TABLE I 
Matrix of Allowed (6,, ¥,,) Combinations Encoded in Binary Logic 


From Table I it is seen, for example, that (B35 ¥3) is chosen when 

0011 = (y > 0, y' > 0, eZ 0, e* € 0) occurs or 

1100 = (y < 0, yt <0, e > 0, e* 0) occurs and that 

(65, ¥1) is not possible. 
In general, it is noted from this encoding that a binary number and its com- 
plement are identical (e.g., 0110 = 1001). Only eight of the sixteen Ona) 
conbinations are allowed. 
The design of a binary logic switching circuit corresponding to Table I 

is straightforward. The signs a the four variables y, y', e and e! are 
sensed and on the basis of the 2+ binary decisions ‘the allowed (6,, yy) com- 
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binations are connected as required. Such a switching circuit has been 
constructed and used in conjunction with an analog computer for simulation 
studies of this type of nonlinear control system. The only essential compo-= 
nents required are relays for switching and zero-coincidence detectors 


ae kes selectors) for sign sensing (i.e., to "read in" the binary logic of 
Table I) 


EXPERIMENTAL STUDIES 


In the original analytical derivation of Eq. (1) it was not possible to 
determine analytically the actual magnitudes of theA , andy, parameters for 
optimum performance. Further, it was assumed that the switching from one par- 
ameter set to another would be done ideally (i.e., instantaneously). However, 
once Eq. (1) has been synthesized, simulation techniques (mentioned in the 
previous section) can be employed to optimize the parameter values and to 
study the effects of the time delay in switching encountered physically. 

A sample of the experimental results is shown in Fig. 2." The response (out- 
put, y, and error, e) of the nonlinear system is compared with the response of 
the linear system (om =Yn = 0). In both cases the input was a triangle wave 
whose frequency was varied from 0.lapt to 0.5wy in Olu, intervals. 


The following data are pertinent to these results: 
Nonlinear system 


D= 0.6 


it 
ine) 


3 = Fy 
po 


Linear system 


iT 
ine) 


¥3 = --Yo 
-fy = 0.5 Yoru 1 


Bhs 


D = 0.6 
Commenting on these results it is noted: 


1. The time averages of the magnitude of instantaneous error, 


[ 
Biz ( terat ° 


compare as follows 


(el nonlinear © 9+1 lel tinear 


where the average is over the entire length of run. 


*Detailed results will be published later. _ | | 
“7. is the natural frequency (rad./sec.) of the undamped linear member. In 
the simulated system, has been chosen to be unity. Ina steady state 
amplitude response of the linear member to sinusoids wy is also essentially 
the cut-off frequency, W,. 


The tick marks at the bottom of Figs. 2(a) and 2(b) denote when the frequency 
was changed. Further, it should be noticed that Figs. 2(a) and 2(b) are two 
separate sets of results in that they were not obtained simultaneously. 
Therefore, there is not exact synchronism in the two time scales. 


fd 


FIG.2(a).-NONLINEAR System Response 


eV 


2. In order to obtain the same improvement in jej by utilizing linear means 
alone, the natural frequency of the system should be increased by the 
same factor aslelis reduced, i.e., 10 times. 


In second order system this requires an increase of loop gain of approxi.- 
mately 102 = 100. 


PRELIMINARY EXTENSIONS TO HIGHER ORDER SYSTEMS 


From a practical viewpoint limitations in applicability of this type con- 
trol do not stem from the feedback switching circuit, but rather from assuming 
that the linear member is second order. In many cases a more realistic 
approach is to consider the linear member to be of higher order but still pre- 
dominantly of second order. A third order system of this nature is currently 
under investigation and promising results are being obtained. 


NONLINEAR COMPENSATION OF AN AIRCRAFT 
INSTRUMENT SERVOMECHANISM BY ANALOG SIMULATION* 


D. Lebell 
The Ramo-Wooldridge Corporation 
Los Angeles, California 


Summary oe Poon of certain analog computer techniques with the direct 
synthesis philosophy* applied to a particularly appropriate servo configura- 
tion’ results in a method for direct specification of compensator design. 

The method centers about a calibration technique which is simple and easily 
applicable to nonlinear control systems of considerable generality. Mathemat- 
ical representation of the physical system is maintained at the realistic 
level typical of analog simulation. 


A detailed description of the method is presented by means of its appli- 
cation to a practical example -- nonlinear compensation of an instrument 
servomechanism. Equations governing the behavior of the unalterable compo- 
nents were written and substantiated experimentally. These equations uniquely 
determine the form of the ideal compensation for the designated system per- 
formance and system configuration. The actual unalterable elements are then 
employed as an analog simulator to "compute" the specified compensator charac- 
teristics. A compensator possessing these characteristic is assembled and 
over-all system performance obtained. An electronic analog computer is 
employed? as a computer simulator for this work. ,; 


The majority of servo compensators designed in the past has been of the 
linear type. It is well known that this partiality is not due to superiority 
of the linear over the nonlinear kind but rather that linear designs have been 
much simpler to perform and nonlinear components hard to come by. The advent 
of analog computers, the extensive progress made in components development, 
plus the increasing need for "getting the most" out of available power trans- 
ducers. combine to encourage removal of the highly restrictive linearity 
condition. 


Nonlinear design techniques available to the engineer have frequently 
suffered from mathematical complexity or computational drudgery on the one 
hand, or unrealistic simplification on the other. Recently, however, consid-~ 
erable progress has been made in the region between these extremes. This 
paper attempts to further such progress and its dissemination. 


An elementary control system representation is shown in Fig. 1. Typi- 
cally, the designer starts with an element G whose characteristics are fixed 
by a combination of design specifications, commercial availability and local 
tradition. His job is to design compensation (or equalization) in order that 
over-all system performance [response of C(t) to R(t)] be satisfactory. For 
the instrument servo and for many other applications it is desired that the 


“This paper was presented at WESCON, the West Coast Conference sponsored by 
the IRE - West Coast Electronics Manufacturer's Association, in August, 1955. 
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Rit B (t) Cit) 


Fig. 1 - Elementary control system. 


output and input be identical over a range of inputs. For the output and 
input to be equal in this configuration, it is necessary for the H component 
to possess the inverse function characteristics of the G component. For exam- 
ple if C = G(B) = B?, it is necessary that B = H(R) = VR, or if C = /Bdt, it 
is necessary that B = dR/dt. The extent to which this condition H = Gu+ can 
be achieved is the extent to which we can approach the ideal system possessing 
no discrepancy between R(t) and C(t). 


DIRECT SYNTHESIS 


The design or synthesis features are exemplified in Fig. 2. Note that 
R(t) = C(t); G, and consequently B(t) are known. Therefore, the input-output 
characteristic of H is uniquely determined at least for any particular R(t), 
and this must be so for any system of any complexity provided only one unknown 
H exists at a time. 


R(t) 


Fig. 2 - Elementary control system with R(t) = C(t) condition imposed. 


Our purpose now is to convert the previous statement into a general pro- 
cedure for obtaining the actual input-output or transference characteristic of 
H. Figure 2 can be mechanized in a computer set-up as indicated in Fig. 3 
where the H position is taken by a recording device (cathode ray oscilloscope) 
whose vertical and horizontal excursions are controlled as indicated. The 
pattern traced on the scope is the solution which the computer set-up yields, 
namely the desired input-output characteristic of H. G™~ may be assembled 
from computer components in accordance with the governing equations. If real- 
time simulation techniques are employed, G71 may be the actual component. 


TECHNIQUES FOR REVERSING SIGNAL DIRECTION THROUGH G 
Some comments are certainly in order with respect to_the practical prob~ 
lems of reversing the direction of signal flow through @-l, If an analog 
computer is used to represent G it is not difficult to invert the transfer- 


ll 


ence. One simply rewrites the differential equation so that the driving 
function and the dependent variable exchange roles. However, if the actual 
physical component G is used (which has the advantage of more realistic system 
representation), its mechanization may present a difficulty. 


One approach is to build an auxiliary high-gain servo (it is convenient 
to use analog computer components here) which compares the actual output of G 
with the desired output, causing the amplified discrepancy to actuate B and 
incidentally the cro. This approach is quite satisfactory within the stabil- 
ity and accuracy limitations of this closed loop around G. 


A second approach has the advantage of increased simplicity (Fig. h). 
A signal is applied at what is normally the input to G [i.e., B(t)]. This 


Fig. 3 - Mechanization of Fig. 2 configuration for direct synthe- 
sis cf compensator (H). 


R(t) = C(t) 


Fig. i - Revision of Fig. 3 arrangement to permit normal signal- 
flow direction through (G). 


signal is chosen so that the output of G conforms to the desired range of R. 
For the moment, G is serving as an analog of itself and also as a signal gen- 
erator for R(t) which is fed to the cro. 


INTERPRETATION OF CRO PATTERN 
What are the typical cro patterns one might expect to see? Suppose H is 


a straight line of slope m. This indicates that H physically should be an 
amplifier of gain m (modified by scale factors as necessary). 
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Suppose H is a curved line which stays substantially constant over the 
signal range. This indicates that a nonlinear transference having this con- 
tour as its specific characteristic is required for H. A compensator with 
this characteristic may be built from nonlinear components such as thyrite, 
germanium, etc. H may display infinite slope in some regions. . This situation 
frequently occurs as a result of saturation or static friction inG. It con- 
stitutes G's legitimate claim on H if the R(t) = C(t) condition is to obtain. 
If the region is significant it may be approximated by a lecal high-gain 
region in the H compensator. 


Alternatively, the cro pattern may appear elliptical (on application of a 
sinusoidal excitation). This is the familiar lissajou figure indicating the 
need for a linear dynamic characteristic for H. 


Finally, a combination of the foregoing may appear, in which case a non- 
linear dynamic characteristic is required. For proper interpretation of the 
cro pattern, however, it must be possible to separate the nonlinear from the 
dynamic characteristic. One must decompose the H component into subcomponents 
each of which is linear-dynamic or nonlinear and nondynamic. The correct | 
order of decomposition must be observed if the inverse nature of G and H is to 
exist. The order in which the H subcomponents are arranged must be inverse to 
that of G. Strictly speaking, the method is only applicable if no more than 
one subcomponent is unknown. By judicious choice of signals, however, it may 
be feasible to apply the method successively to each part of H. For example, 
a high-frequency low-amplitude R(t) will permit investigation of the linear- 
dynamic block; a low-frequency high-amplitude signal will isolate the nondy- 
namic nonlinearity. General knowledge of the nature of G will aid inestimably 
in adapting the method to new situations. Where decomposition of H yields 
more than two unknown subcomponents the method does not appear to be gener- 
ally practical. 


THE UNALTERABLE ELEMENT G 


Conforming to recommendations of the preceding section, preliminary study 
of unalterable element G is helpful in establishing the form of H to be anti- 
cipated. For the particular magnetic amplifier-motor-gear train combination 
employed the torque equations are: 


C25 4)mq) BE IC a): 
where 


-1 for T< 0 

i for Tt. 0 

developed torque 

effective moment of inertia of motor and load 


am a 
Te | ed | 


per-unit motor acceleration 

amplitude of two-phase 00 cycle voltage 
constant of proportionality 
1.0 for suitable choice of units. 


wa 
Ut | 


U 
oO 


* _ 24 2 
valid for at 
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Note that developed torque is assumed proportional to the praduct of slip 
and control voltage and that representation of slip for bidirectional motor 
operation requires introduction of the sign operator q. Also, the assumption 
which more or less defines an instrument servo is invoked, namely, that fric- 
tion torque is negligible in comparison with inertia torque. Finally, the 
magnetic amplifier time delay is negligible. 


CONFIGURATION FOR H 


Equation (1) leads to 


JC = B(l - qc) = Ba (= -C) (2) 
or 
= JC - 2d J 75, 
. a(1/q-6¢) dtq ea 
since 
q=1/q 


In what manner must B be generated from R for the condition (R = C)? 
Replacing C with R in Eq. (2) results in 


Susiedlyypenet * 
fey are aed In(q - R) (3) 


wnich defines the form of H. This relationship is further depicted. in the 
diagram in Fig. 5. i 


q GENERATOR 


H 


Fig. 5 - Representation of compensator (H) from Eq. (3). 
as A dq 
valid for aE 6) 
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INPUT 
SIGNAL 


G 


(MAGNETIC 
AMPLIFIER, MOTOR 
GEARING, ETC.) 


ACHOMETER 


N | DETERMINATION 


Fig. 6 - Experimental set-up for cro display of nonlinearity. 


Fig. 7 - Obtained 
from set- 
up of Fig. 
6 is cro 
pattern. 


Ne 


EXPERIMENTAL SET-UP 


Figure 6 illustrates the arrangement for producing the nonlinear charac- 
teristic of the compensator. The motor tachometer serves conveniently as one 
of the H differentiators. The other is temporarily an integrator since the 
direction of signal flow is reversed. It is constructed of analog computer 
elements. 


Finally, note the deletion of the q generator. This last simplification 
produces a piecewise-constant error (known) in the cro pattern which is easily 
corrected. Theory predicts that the resultant cro pattern should appear 
piecewise logarithmic which it does (Fig. 7). This pattern was found to be 
essentially invariant over the operating range of input signals. 
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INTERPRETATION OF NONLINEAR CRO PATTERN 


The multivalued nature of the curve while aesthetically pleasing requires 
some interpretation. The tachometer output is a 00 cycle voltage whose peak 
magnitude (envelope) is proportional to speed. For speed reversal, the tachc- 
meter voltage reverses phase with respect to line voltage, a phenomenon which 
is lost in a high persistence cro. From the expected system behavior, how- 
ever, one can easily select the appropriate half of the envelope. The 
function is now reduced to a mere double-valued curve. This last ambiguity 
was to be expected since, as mentioned earlier, the abscissa of the cro pat- 
tern is R rather than the q - R which would have produced a single-valued 
function. Replotting the curve as a function of q - R leads to the N71 char- 
acteristic eventually constructed from computer components. Note that we may 
select an arbitrary horizontal datum for the two branches of the curve sinee 
the output of N-L is to be differentiated. This was done to permit the _ 
smoothest approximation to the function. In so doing, two biased diodes suf- | 
ficed to approximate the curve in the final set-up made on the computer. 


SYSTEM CONFIGURATION AND PERFORMANCE 


The simple configuration in Fig. 1 suffers from the absence of feedback 
and feed-forward. The configuration in Fig. 8 was chosen to desensitize sys- 
tem performance to component variations. In so doing, it was possible to 
employ low tolerance components (10 - 25%) and yet obtain the performance 
indicated in Fig. 9. The data compare performance of the nonlinear, dynamic 
compensator (Fig. 5) with that of a conventional, lead-lag compensator. The 
comparison is based on an average error (divided by input) taken over a range 
of input amplitudes at each frequency for which a data point is indicated. 


Several items may be noted with respect to Fig. 9. The sharply rising 
error (as frequency increases) is the typical "cut-off" phenomenon. An 
improvement in high frequency performance is attained through a combination of 
two factors: : 


1. The nonlinear characteristic in H compensates in part for the nonlinear 
saturation inherent inG. 


Fig. 8 - Final system configuration employed. 
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Fig. 9 - Comparison of final system performance with 
that of conventional lead-lag compensation 
(direct feedback). 


2. The configuration of Fig. 8 permits a greater "lead" or anticipation in 
H to compensate the "lag" inG. This is due to the stability limita- 
tions inherent in the conventional direct feedback configuration. 


On the other hand, since the lead-lag compensator permitted closure of 
the feedback loop at a higher loop gain than in the other case, it displays 
slightly superior performance in the low frequency region. (Note that H does 
not pass low frequency signals.) It is apparent that a system employing both 
types of compensation will possess the attributes of both. 


It will be noted from the preceding paragraph that as is frequently the 
case, improved performance comes at the price of increased complexity. The 
usual principles of good design demand that a proper balance be struck between 
complexity and performance for each situation. It should be emphasized that 
the techniques illustrated are largely independent of the over-all system com- 
plexity. They have been tested successfully on other systems of practical 
interest. 


CONCLUSION 
Restrictions on the method described are: 


1. The unknown compensator (EK) must be representable (at worst) by two 
unknown subcomponents -- nonlinear-nondynamic and linear-dynamic. 
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2. The cro patterns which result must be capable of synthesis from physi- 


cally realizable components. Otherwise, H is incorrectly located, or 
the performance criterion cannot be satisfied. 


3. A particular cro pattern is valid for the particular R(t) excitation 


only. It must remain approximately invariant over the range of antici- 
pated signals if it is to be synthesized from nonlinear-nondynanii.c 
elements. 


If these conditions can be met, a straightforward calibration method of 


considerable generality is available wherein the unalterable elements of a 
system prescribe directly the compensator required to meet the performance 
specificaticn. 
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A-STEADY STATE APPROACH TO THE THEORY OF SATURABLE SERVO SYSTEMS 


J.C. Lozier 
Bell Telephone Laboratories 
Whippany, New Jersey 


Summary -- It is shown that simple saturable servo systems can have two modes 
of response to a given input signal -~ one a linear mode as predicted by lin- 
ear feedback theory, and the other a saturated mode predictable from a large 
signal analysis presented here. This dual mode of response theory provides a 
reasonable explanation for premature saturation, hysteresis in the input- 
output characteristics, jump resonance and similar anomalous effects that so 
often degrade the steady state performance of saturable servo systems. Fur- 
thermore, systems which exhibit a large degree of anomalous steady state 
behavior can be expected to exhibit a correspondingly poor transient response 
under large signal conditions such as obtained when the system is coming out 
of saturation. 


Demonstration of the existence of this dual mode of response and predic- 
tion of the range of signal amplitude and frequencies for which it exists are 
both built around the "saturated" transfer characteristics of the control 
loop. This analysis is quantitatively useful only in certain simple cases 
where the necessary saturated transfer characteristics can readily be found. 
Its chief value lies in the insight it provides to the cause of anomalous 
behavior. In this respect it is particularly useful to the designer because 
it relates this behavior to the loop gain and phase characteristics of the 
saturated transfer characteristics. 


The analysis is applied here to two simple examples, and the results are 
verified by an analog computer study on simulations of these systems. 


The type of saturation discussed here is the familiar power saturation 
that occurs in the forward path of a linear feedback control system and marks 
the end of its linear range of operation. In a system where size, weight, 
efficiency, etc., are important considerations it is often necessary to com- 
promise on the power handling capacity to the point where operation under 
power saturated conditions is to be expected a good part of the time. Such 
systems will be referred to here as "saturable" servo systems. 


The performance of these systems obviously depends largely on the per- 
formance in the neighborhood of saturation. Unfortunately, servo systems 
often exhibit a large degree of anomalous behavior in this region. This anom- 
alous behavior can affect the steady state performance by causing premature 
saturation, hysteresis, "jump resonance"’*“ and similar double valued system 
response characteristics, Enhancement of the amplitude and lengthening of the 
period of the overshoot, 2 reduced stability margins in the presence of noise 
and even low frequency oscillations are transient manifestations of such anom- 
alous behavior. 


All the anomalous effects noted above are large signal phenomena. The 
problem of how a saturable system will perform in the presence of large sig- 
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nals is a nonlinear one. A linear feedback analysis of the system is helpful 
to the extent that a satisfactory small signal performance is a prerequisite 
to satisfactory large signal performance. However, linear concepts are often 
misleading for handling large signal performance. 


Fortunately a form of steady state analysis of these quasi-linear systems 
can generally be made by fixing the amplitude of the signal at each frequency. 
This is the basis of Kochenburger's analysis of contactor servos,+ and it is 
the basis of this analysis. This approach will be used to derive what is 
termed here the "saturated" control loop transfer characteristics for these 
saturable systems. With this saturated transfer characteristic, the possibil- 
ity of either a linear or a saturated response for the same input signal can 
be demonstrated. The conditions for the existence of this dual response 
behavior will be shown to depend on the gain and phase of the saturated trans- 
fer characteristics. 


The technique is illustrated with two simple examples and the computed 
results are verified by the use of an analog computer. 


DESCRIPTION 


In the course of the analysis, two feedback control systems will be stud- 
ied and compared. They will be referred to here as the "reference" system and 
the "modified" system. The following assumptions and restrictions concern 
these systems. 


Assumptions and Restrictions 
1. Both will be simple single loop systems of the type shown in Fig. l. 


2. In Fig. 1 the box labeled "servo amplifier and equalizer" will com- 
prise the control elements. It will be assumed that the transfer function of 
this box, Go3> will be independent of the amplitude of the input signal. 


3. The only difference between the reference system and the modified 
system is that Go3 is just a constant factor K in the case of the reference 
system, whereas G23 will have a specified amplitude and phase vs. frequency 
characteristic, as well as a gain factor, in the case of the modified system. 


h. The only nonlinearity in these systems will be in the box labeled 
"limiter." This limiter will have an output amplitude proportional to its 
input amplitude up to the limiting point, as shown in Fig. 2. The limiter is 
assumed to have no amplitude or phase vs. frequency characteristics, so that 
its transfer function G3), is just a constant gain factor below limiting and a 
variable gain factor above. 


5. The box labeled "servo motor and load" will be assumed to have a lin- 
ear input-output characteristic for all input amplitudes passed by the lim- 
ter. In other words, the saturation takes place in the limiter just ahead of 
the motor. 


6. The steady state transfer characteristics of the motor and load, G)\5, 
is assumed to be that as shown in Fig. 3. 
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Note: The shape of this transfer function is fairly typical of such motors. 
The constant velocity and constant acceleration asymptotes which 
bound the output displacement vs. frequency characteristic over much 
of its useful frequency band follow directly from the assumption that 
the velocity and the torque the motor can develop are proportional to 
the amplitude of the input voltage. A parasitic asymptote such as 
the one assumed here is also typical of such motors. It was added 
for completeness and has little bearing “on the main arguments given 
here. The associated phase vs. frequency characteristic assumed here 
is simply the minimal” phase that must accompany the assumed ampli- 
tude vs. frequency characteristic. 


7. It will be further assumed that the shape of the displacement vs. 
frequency, and of the phase vs. frequency characteristics of the motor are the 
same whether the motor is driven by a sine wave signal or by a square wave of 
the same frequency. 


Note: This assumption ignores the harmonics of the fundamental component 
that may be present in either the input or output signals. Such 
assumptions are customarily justified in computations of this type on 
the basis that feedback control systems are low pass devices. One 
exception to this will be made below. : 


8. It is assumed that neither the relative shape of the displacement vs. 
frequency characteristic nor the actual phase vs. frequency characteristic of 
the limiter plus the motor is changed by saturation. In other words, satura- 
tion in the limiter and motor combination is assumed to just change the flat 
gain factor in the limiter transfer function. 


These assumptions have been chosen to make it easy to derive the steady 
state saturable transfer characteristics for the limiter plus motor and load. 


Derivation of the Saturated Transfer Characteristics 
The maximum output displacement vs. frequency characteristics of the lim- 
iter plus motor and load characteristics may readily be derived from the 
linear transfer characteristics shown in Fig. 3, using the above assumptions. 
In the linear case, for an input voltage to the limiter of the form 
e(t) = E cos wt 

the output displacement of the motor and load is 

Go(t) = [®o| cos (wt + 9) 
and Fig. 3 shows curves of 20 logig Pc and % as functions of frequency. 


In saturation the output contains a fundamental component 


l(8o)ml cos (wt + Z) 


“Bode, H. W.: "Network Analysis and Feedback Amplifier Design." D. Van 
Nostrand Company. 
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Female 
with | (85) ml constant at any given frequency for any input voltage E exceeding 


the saturation value E,, Figure l shows a plot of 20 ogy Len! as a func- 


tion of frequency, where K is a constant. This maximum output displacement 
vs. frequency characteristic has the same relative shape as the linear charac- 
teristic due to assumptions h, 5, 6 and 7 above. The saturated phase 
characteristic is identical to the linear phase characteristics again because 
of assumptions h, 5, 6 and 7. 


Based on these simplifying assumptions and restrictions, the desired 
steady state saturated transfér function for limiter and the motor plus load 
can be expressed as 


(85) m(w) = G30(w)E.. 61545) 


Here (@5)m is the maximum output displacement vs. frequency, and E, is the 
input to the limiter that will just produce saturation. One major difference 
between the saturated transfer function given in Eq. (1.1) and its linear 
counterpart is that this saturated transfer function only holds for one single 
input frequency” at a time. This follows from the fact that the superposition 
theorem used in linear analysis to handle multiple input signals breaks down 
in this case. 


In many saturable transfer elements, both the shape of the amplitude vs. 
frequency characteristic, and the phase vs. frequency characteristics will 
change with the amplitude of the saturating signal. In these cases, of 
course, the saturated transfer characteristics would be different from the 


"By assumption 7, repetitive square waves or triangular waves are included. 


23 


linear characteristics. In fact, there may actually be a whole family of such 
saturated transfer characteristics, depending on the precise amplitude of the 
input signals. Such factors would complicate application of the steady state 
analysis presented below. 


STEADY STATE ANALYSIS 


The first part of this section will be devoted to showing the basic tech- 
niques used here. The effects of output saturation in producing hysteresis in 
the steady state input-output characteristics will be demonstrated, and the 
equations governing the relationship between hysteresis and the saturated loop 
transfer characteristics will be derived. 


In the second part of this section, the technique will be applied to two 
illustrative examples. 


The Basic Technique 


The technique used here is based on the relationship that must exist 
between the three signals associated with the linear summing point in the 
servo shown in Fig. 1. These signals are the loop input signal x;(t), the 
loop feedback signal x,(t) and the difference or error signal xp(t). Since x» 
is the difference between xj and xp, the equation governing this relationship 
is simply 


X_(t) = x(t) - x(t). (1) 


Now for a single frequency input, x;(t) = A; cos (wt), x,(t) and x,(t) are 
also sinusoids of the same frequency, ignoring harmonics generated by the sat- 
uration in the output. Hence, we can transform Eq. (1), which is good at any 
instant of time, into the vector equation, 


Xo(w) = X(w) - X,(w) (2) 
which holds good at any frequency w. 
According to Eq. (2), these three vectors must form a triangle at any 


frequency w. The triangle that is most useful for illustrative purposes is 
that shown in Fig. 5. In drawing the triangle in this way, X, rather than X, 


is the reference vector. This is convenient because X, and X,are related by 
the loop transfer function G, so that 


Xo(w) = G(w) Xo(w). (3) 


Hence, with Xe as the reference vector, -Xp can be readily determined from G. 
Once Xe and -Xo are drawn as shown, then Xj must be the vector connecting =X, 
and Xe, Incidentally, this method of representing Xe and -Xo is related to 
the one used in drawing the Nyquist diagram of a system. It should be noted 
that the angle between X, and =X, is 180° - /G and therefore represents the 
phase margin of the control loop transmission at the frequency in question. 
The symbol @ will be used for this phase margin, so we have 


@ = 180° - /G. (1) 


Introduction of Output Saturation 


The main interest here centers on what happens as Xp saturates. From the 
assumptions made earlier, it follows that Xo reaches a maximum value (Xo)m 
when = reaches the value (Xe)s. Xe is free to increase beyond the saturation 
value (Xe)s, but -(Xo)m and /Gs, the angle between Xe and -(Xo)m, do not 
change. Hence, as the input Xj is increased beyond (X;),5, the level that pro- 
duces saturation, the magnitude of X,_ must change in order to complete the 
triangle of X,, -(Xo)m and X;. Figure 6 shows this phenomenon graphically. 

In Case (a) the input is below saturation; in Case (b) the input is at satura- 
tion levels; in Case (c) the input is above saturation. For comparison 
purposes, Fig. 6 also shows what would happen in the equivalent linear case, 
where no saturation occurs. As Case (c) shows, saturation prevents -Xg from 
increasing beyond - (Xo)m, hence Xg, the error, must increase as X; increases 
to maintain the vector relationship between -(X,),, X; and X,. This example 
in Fig. 6 shows what happens in a well behaved case of saturation, where the 
relationship between X; and -X, is single valued, and all three vectors can be 
conclusively determined. 


The example shown in Fig. 7 is one where the large signal behavior is not 
so well behaved. We see here the same three input conditions that were used 
in Fig. 6, and the equivalent linear case is again drawn for reference pur- 
poses. In Case (a), where the input is small, the system is effectively 
linear as before. In Case (c), we again have no trouble determining the vec- 
tor relationship. However, in Case (b), where the input level is just enough 
to produce saturation, two possible angular positions are shown for Xj, and 
two possible magnitudes for X,. 


Conditions for Anomalous Large Signal Behavior 


Further examination of the example in Fig. 7 will show that (X;), is not 
the only input amplitude for which the input-output characteristics of the 
loop has double valued possibilities. Another case is shown in Fig. 8. Here 
the loop characteristics are again the same as those for Fig. 7, although the 
scale factor of the drawing has been expanded for clarity. The input ampli- 
tude shown here is intermediate between cases (a) and (b) shown in Fig. 7. 
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Seip eo 


The precise input amplitude used in Fig. 8 is designated (Xz) and will 
be called the "critical" value. The analysis below will, show why this criti- 
cal value is of particular interest. 


As Fig. 8 shows (XS Jel is the magnitude of the perpendicular from (Xo)m 
to the Xe axis. Hence, we can write 


(X4)o| = (Xo) sin’d. (5) 


The importance of this critical input amplitude stems largely from the 
fact that it is the smallest amplitude of X; that could possibly form-a closed 
triangle with X,, © and Tea as shown. Hence, no input amplitude 
\Xi1 < |(Xi)e] could produce saturation, and the system must be linear for 
these smaller amplitudes. 


However, when [X;| =|(X;),| the system can have two possible responses; 
one is the linear response (Xe and the other is the saturated response 
(Xo)me The linear response is obviously a stable operating condition since 
{(X4)e] < 1(X%:)5].. As for the saturated condition, (Xj)¢~ can form the neces- 
sary triangle with (Xo)m, Xe and @ as shown in Fig. 8. The question of 
whether this triangle represents a stable operating condition depends on 
whether or not the magnitude of X, involved is large enough to sustain the 
output at the saturation level (X_),. Now the magnitude of X, involved is 
readily seen to be 


[Xe] = |(Xom| cos 0. (6) 


The minimum amplitude of X, that can sustain the output at saturation is 
(Xe)5. This triangle, therefore, will represent a stable operating condition 
when 


(ear < [onl cos 0. (7) 
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In the example shown in Fig. 8 this inequality is satisfied, and the sec-~- 
ond triangle with |X| = peas is a possible operating condition. 


Thus when Eq. (7) is satisfied, the system has two possible equilibrium 
conditions for the same input amplitude |X;| = |(X,),|. This represents a 
type of anomalous behavior for which the output response is either the linear 
value corresponding to the input amplitude or the saturated value Conia 


Similar arguments can be presented to show that such anomalous large 
signal behavior must exist for all input amplitude between this critical input 
amplitude |(X;)¢] and |(X;)5], the input amplitude for which saturation is to 
be expected. 


The conditions for anomalous large signal behavior at saturation expres- 
sed by the inequality in Eq. (7) deserve further discussion. It was assumed 
that the system under discussion was linear fcr al] error signals 
((Xe)] < |\(Xe)sl|. Hence, we write 


(Xo)m = G(Xe) se (8) 


When Eq. (8) is combined with Eq. (7), we see that (Xe) s| < |G||(%e)s] cos 6 
or, cancelling (Xe)s, that 


|G| cos 9 > 1. (9) 


This equation expresses the conditions of anomalous steady state behavior 
in terms of the loop transfer characteristics. 


Effect of Saturation on the System Transfer Characteristics 


The steady state operation of a simple linear servo system can generally 
be described by three system transfer functions relating the three variables 
X;, X_ and X,. These three transfer functions are: 


1. The steady state input-output transfer function, relating X, to Xi, 
which can be written for the linear case as 


= eet (10) 


2. The steady state input to error function, relating X, to X;, can be 
written for the linear case as 


x 
££: ar Ste “115 
Xi 1+G 


3. The loop transfer function G relating X, to X_ can be written for the 
linear case as 


—=G. (12) 


The large signal analysis wil] now be used to determine how these three 
linear transfer functions are modified when saturation is present. As the 
previous analysis shows there are two concitions of interest. 


The first condition is the one for which the locp transfer characteris- 
tics are such that 


IG| cos 9 < 1. (13) 


In this case, no anomalous large signal will occur, and the linear Eqs. (10), 
(11) and (12) still hold for the three transfer functions. However, their 


applicability must be restricted to input amplitudes (3 | Z | (Xi) s| because . 
of saturation. 


The second condition is one for which 
Ic] cos@>1. (9) 


In this case, the analysis shows that anomalous behavior will occur for input 
amplitudes above the critical amplitude, \Xi]c, where 


([%le = |(o)m| sin 0. fnck®) 
Hence, there are two subsidiary cases depending on the input amplitude. 


Case a. For |Xs] < |(X;)o], the linear equations still hold for the 
three loop transfer functions, as in the first case above, except that the 
restriction on input amplitudes is determined by (X;)¢ instead of (Xi). 


Case b. For inputs where such anomalous large signal behavior occurs, 
namely for (X;) so that [Xe] < [Xi] < [ (Xi) 5], we have a new set of 
transfer functions. 5 


(1) The input-output transfer characteristics, relating X, to 
;, is now double valued with the output having either the value given by 
the linear Eq.—{10) or the saturated value |(Xo)m|. Hence, we have the 
anomalous situation where 


Xo G 
oe ee ee ee 
or | (14) 
IXol \o)mI * 
[a] | a | 
(2) The steady state input to error relating X, to X, is also 
double valued in this case. It can have the linear value 


b] 


bi 
@ 
KB 


*?he angle between Xo and X; is ambiguous. 
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or, when the output is saturated, 


% 


2 
| xel 3 1 (x4 )el antes: —_ \(xidel . 
Pal | C%a) | \(x4) |? 


(3) The loop transfer function relating ato Xe is unaffected 
by this anomalous behavior. ‘The linear relationship given by Eq. (12) 
still applies for amplitudes of | Xe] < | (Xe) s] - 


(15) 


EFFECT OF SATURATION ON THE STEADY STATE PERFORMANCE 


Using the system transfer characteristics developed in the previous sec~ 
tion, the effect of saturation on the steady state performance can now be 
discussed. 


The steady state input-output transfer function will exhibit anomalous 
large signal behavior unless: 


|G] cos @ < 1. (13) 


This means that the phase margin would have to approach 90° when G » the loop. 
gain, is large, in order to avoid anomalous behavior. Such a restriction is 
too severe for most practical servo systems to meet at all frequencies, when 
other requirements such as bandwidth or speed of response are considered. 
Hence, anomalous large signal behavior is to be expected over some part of the 
frequency range of the system. : 


The type of input-output response that is obtained under anomalous condi- 
tions is shown in Fig. 9. It shows the typical double valued output that 
occurs for input amplitudes between | (X;)oJand |(Xi),]. This double valued 
effect similar to hysteresis and is certainly not desirable. It is apparent 
that the dependable range for linear operation ends at X; = eure 


The effectiveness of the system in utilizing the available power hand- 
ling capacity under these conditions is measured by the ratio of | (Xq)Qfto 


I(X:) 51, or still better by the ratio of |(X4)cl to |(Xo)m]. By rearranging 
Eq. (5) we get 


[(Xael _ 4 


Gane ee 


“The derivation of this expression is straightforward but lengthy. It is 
based on the fact that 


[xef I(Xo)ml_cos @ Cees 
[ (Xie | Xo)m|sin @ ; 


The phase is again ambiguous. A negative sign in front of the radical in Eq. 
(15) is also possible over a small range of values of X;, but it seems 
unlikely from physical considerations. 
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These results indicate, for example, that at a frequency for which the phase 
margin is 30° and the loop gainG > 1, the input-output response of the servo 
is single valued and linear over only half of its potential range. The impor- 
tance of keeping 9 as large as possible over the working frequency range of 


the system is apparent. 


The steady state input to error transfer function of the system also 
exhibits large signal anomalous behavior under the same conditions as the 
{4 put-output transfer functions. Figure 10 gives an example of such behavior. 
A case of particular interest is that for which |X;] = [(%;)~].. The corre- 
sponding magnitude of X, can have either the linear value 


|Xel > I(Xe)el = oa . (17) 
or | NP 
{xe| = |(%)el cot 0. (18) 


It may be shown that the larger the phase margin 8 at the frequency in ques- 
tion, the smaller the difference that will exist between the two values for 
|X], as given by Eq. (17) and (18). Similarly, it can be shown that for the 
same phase margin 8, this difference increases as JG| increases. This can be 
seen from a comparison of Eqs. (16) and (18). Again the importance of keeping 
@ large when|G| is large is emphasized. 


APPLICATION OF STEADY STATE ANALYSIS © 


In this section we will analyze and compare the steady state performance 
of two systems which will be referred to as the reference and the modified 
systems. -These examples have been chosen in order to demonstrate the conclu- 
sion that in saturable servo systems the loop phase margin is an important 
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factor in the performance, even at frequencies where the loop gain [G| » 1. 
The reference servo has what may be termed a fairly satisfactory loop phase 
margin under high loop gain conditions, and the modified servo has a rela- 
tively poor one. In other respects the two systems are as much alike as 
possible. To make the analysis simple, the seven assumptions made previously 
will apply to both systems. Both will be simple single loop systems as shown 
in: Figs 23 


The first step will be to present the loop transfer functions for these 
two systems. In the case of the reference servo, it is assumed that the con- 
trol element, representing the servo amplifier and equalizer, has a constant 
nonfrequency sensitive transfer function. Hence the linear loop transfer 
characteristics are the same as the characteristics shown in Fig. 3 for the 
limiter plus the motor and load, except for a constant gain factor. The 
resultant loop transfer characteristics for the reference servo are shown in 
Fig. 11, The absolute gain around the ee in this case was chosen so that 
the gain margin at phase crossover (6 = 0) is about 17 db. 
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The linear transfer characteristics for the modified servo were obtained 
by adding an integrating network of the type shown in Fig. 11 to the reference 
servo, and then readjusting the flat loop gain until the gain margin at phase 
crossover is again about 17 db. The resultant loop transfer characteristics 
are also shown in Fig. 11, The constants of the integrating network were 
chosen so as to increase the relative loop gain at frequencies below 0.0078 by . 
2 db with the effect tapering off to zero at frequencies above 0.25 radians 
per second, as shown.. 


Since both the reference and the modified servos use the same limiter 


plus motor and load, the saturated output amplitude vs. frequency characteris- 
tie given in Fig. holds for both systems. 
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Input-Output Performance 


Both of these systems were found to satisfy the conditions for anomalous 
large signal behavior for all frequencies of interest below approximately 0.6 
radians per second. Hence, to determine the effect of saturation in the 


X 
input-output behavior for these systems, we computed can over this fre- 
o/m 


quency range below 0.6 radians per second. Equation (16), which shows that 


: . 
oe = sin @, was used in this case. The necessary values of |(Xo)m| used 
o/m 


in the computations are given in Fig. for both systems. The phase margin @ 
in each case was derived from the linear transfer characteristics in Mig. ll. 
This convenient procedure is allowable on the assumption that 6 does not 
change with saturation in these two systems. The results of the computations 
are both given in Fig. 12. The curve for the reference system shows that . 


(Xa) | 


1G: ) | is very nearly unity except for frequencies in the neighborhood of 
o/m 


w = 0.6 radians per second. This means that the range of input amplitudes 
exhibiting anomalous behaviors is negligibly small, and the system will be 
linear for input amplitudes practically up to |(Xo)ml. On the other hand, the 
curve for the modified servo shows that this ratio drops as low as -7 db at 

w = 0.03 radians per second, and does not approach the ideal value of unity 
until w drops below 0.002 radians per second. This means that the anomalous 
behavior exists for over half of the nominally linear range of input ampli- 
tudes for frequencies around w = 0.03 radians per second. Hence, the input- 
output characteristics ceases to be linear when the output is still 7 db below 
its saturation value at this frequency. 


Compared with the reference servo, this is poor behavior. This poor 
behavior must be due to the reduction in loop phase margin, resulting from the 
introduction of the integrating network. The only other difference between 
the two systems is in loop gain at the low frequencies, and gain does not 
affect this ratio, as Eq. (16) shows. 

| 

| 


mes 
Figure 12 also, shows a plot of the ratio ion computer on the assump- 
m 
tion that the modified servo behaves linearly for input amplitudes below 
saturation. This would be the performance of the system if there were no anom- 
alous large signal behavior. 


A similar comparison of Oa) sh ana Lael for the reference servo 
| (Xo) m| | (Xo) ml 
showed them to be practically identical, which is not surprising in this case 
where the anomalous effect is very small. 


The name "jump resonance" is given to a discontinuity that may occur when 
the input vs. output characteristic of a servo is measured by a sinusoidal 
input signal of constant amplitude that sweeps back and forth over the fre- 
quency range of interest. The phenomenon is distinguished by a hysteresis 
effect in the output caused by a sudden jump from the linear to the saturated 
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output response at a relatively high frequency as the sweep frequency is 
increasing, and a continuance of the output at the saturation level during the 
return sweep until, at some lower frequency, the output falls back to the lin- 
ear value. This phenomenon occurs only when the input amplitude falls within 
certain limits. Such a hysteresis effect in the input-output characteristic 


\@el , 
(eal? 
small, because a signal amplitude that equals the linear saturation value 
[(x4)5| at one frequency would not fall below | (X)_] until a lower frequency 


was reached on the reverse frequency sweep. Under such conditions, a hystere- 
sis of the type ‘described is quite reasonable. 


with frequency could occur in the case of the modified servo where —————— 


Input-Output Measurements 


YU 
The steady state results on these two systems can be checked, so an ana- 
log computer was set up to simulate these two systems. Then the output 
amplitude was measured as a function of the input amplitude. The input fre- 
quency used was 0.0316 radians per second. The reference servo showed no 
measurable anomalous behavior at this frequency. The modified servo gave the 
results shown in Fig. rsa) This measured input-output characteristic is dual 
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valued, as the theory predicts. Furthermore, the amplitude range over which 
anomalous behavior exists was measured as approximately 7.8 db. This checks 
closely with the theoretical results which give 7.2 db. 


Note: It was found necessary when checking these steady state conclusions 
on an analog computer, to make allowance for the fact that the satu- 
rated output displacement of the motor is a sawtooth rather than a 
sine wave at relatively low frequencies. This is because only the 
maximum velocity limitation is important at low frequencies. Under 
these conditions the feedback voltage is also a sawtooth. This makes 
it difficult to explore the steady state performance of the system in 
the immediate neighborhood. of saturation with sinusoidal input sig- 
nals, because the difference or error signal between, a sinusoidal 
input signal and a sawtooth feedback signal of approximately equal 
amplitude has a large third harmonic term. As a result, under some 
conditions the error can change sign three times in a half cycle of 
the fundamental. To overcome this difficulty, a sawtooth rather than 
a sinusoidal input signal was used. It was not expected that this 
would substantially change the results since the amplitude of the 
third harmonic in a sawtooth is 19 db below the fundamental. This 
expectation is borne out by the fact that the results of the computa-~ 
tions assuming sinusoidal inputs check very closely with the analog 
computer measurements using the sawtooth input signals. ; 


Input-Error Performance 


The steady state analysis shows that anomalous behavior in the input- 
output characteristics of a saturable system is accompanied by a double valued 
input-error characteristic. These results are more difficult to analyze 
because the error signal does not stay fixed under saturated conditions as 
does the output. This is shown by Eq. (15). 


Therefore, for comparison purposes, we have computed the error to input 
signal ratio for the modified system at the critical input amplitude ((Xaels 
under the assumption that the output amplitude has jumped to its saturated 
value, and under linear conditions. The equation for the error under satu- 
rated conditions is given by: 


= cot 8. (19) 
[CG] 
Ixel Se 
Figure 1 shows this ratio * Tae for the modified servo system. gure 
i/e 
ly also shows the linear case where 
1 faa ale 


i i ies above 0.6 radi- 
The modified system becomes effectively linear at frequenci 26. 
ans per second, where anomalous behavior ceases, and the linear equations were 
used to extend the analysis up to 2 radians per second. The reference servo 
gave results similar to the linear case in Fig. ll. 
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As a comparison of the curves in Fig. 1) will show, the error to input 
ratio is greatly increased by anomalous behavior in the case of the modified 
servo, with a peak greater than unity occurring roughly at w = 0.03 radians 
per second. This corresponds to the dip in effective load carrying capacity 
that occurs at the same frequency, and for the same reason. The phase margin 
has a minimum at this frequency. The undesirability of an enhanced error 
characteristic of the type shown for the modified servo system is self 
evident. 


Input-Error Measurements 


Measurements of this hysteresis effect on the input vs. error character- 
istics were also made on the analog computer for the modified servo at 
w ® 0.0316 radians per second. These results are shown in Fig. 15. The com 
puted jump in error at the critical input amplitude is 6 db and the measured 
value is 5.5 db. 


TRANSIENT CONSIDERATIONS. 


The large signal analysis presented above is restricted to single fre- 
quency signal conditions and so is not applicable to transient problems. The 
linear theory is not applicable either when the transients are large enough to 
produce saturation in the system. However, some useful information can be 
obtained from studying specific cases. 


In a saturable servo system a large input step function will initially 
reduce the gain around the loop through compression in the limiter. The gain 
will recover in time as the output approaches the new equilibrium position 
required, and the fed-back signal restores the balance at the input. Hence, 
the gain is changing under these transient conditions. If the system is one - 
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like the reference system analyzed above, in which the anomalous steady state 
behavior is negligible, the form of the transient recovery can generally be 
determined from a knowledge of the maximum velocity and the linear transient 
response to a step function just large enough to produce saturation.* When 
the system is like the modified system analyzed above, the anomalous steady 
state behavior is severe, and it is not likely that the transient could be 
represented by a simple damped exponential. However the transient recovery 
from a large step function is likely to be seriously degraded. Certainly the 
transient behavior would be less damped, and frequency of the oscillating com- 
ponent would be lower if the loop gain were artificially held fixed at a 
substantially reduced value. 


Transient Measurements 
SS a Se 


The conclusion that a system exhibiting severe anomalous steady state 
behavior will also exhibit a degraded transient behavior in recovering from 
saturation is borne out by the transient response measurements made on an ana- 
log computer set up to simulate the reference and modified servo systems. 
Figure 16 shows the results for the two systems for a step function input. 
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*16 satisfy this requirement for negligible anomalous steady state behavior 
-for |G| large, the minimum phase margin should be 60 
Li (X4)e| = 0-87 |(X—)m|]- Then the transient behavior at reduced gains would 


not be under damped. 3 
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The step input was 5 times as large as the level that would just saturate the 

system. The reference servo showed an overshoot of approximately 1/10 of the 

linear range of operation, and the transient died out is less then 6 seconds. 

The modified servo, on the other hand, had a transient overshoot 2 4 times the 
linear range of operation, required aly seconds to return to the balance point 

and then showed a second and third overshoot. 


Stability vs. Noise 


The phenomenon of a servo system becoming less stable in the presence of 
noise is readily explainable if the system is a saturable one in which the 
phase margin is substantially less at frequencies below gain cut-off. Then 
noise of sufficient amplitude to produce limiting would reduce the loop gain 
to desired signals and thus effectively reduce the phase margin of the system 
at gain crossover. Furthermore, if the noise were of much higher frequency 
than the signals, the amplitude of the noise required to produce velocity or 
acceleration limiting would be correspondingly less than that of signals. 

This follows from the fact that the load carrying capacity of the system falls 
off with frequency. 


DESIGN IMPLICATIONS 


The emphasis here is on understanding this anomalous behavior in satura- 
ble systems rather than on their design. However, it should be recognized 
that to eliminate anomalous behavior by restricting designs to ones with suf- 
ficiently large phase margins at high loop gains is not realistic. For other 
reasons it may not be desirable or even possible to do this. It is more real- 
istic to treat these gain and phase indications of anomalous behavior as a 
warning of when and where these troubles will arise. Other techniques for 
- minimizing this anomalous behavior can then be developed where necessary to 
fit particular cases. For example, in the modified servo analyzed above, it 
should be possible to greatly minimize the anomalous effects by the use of a 
limiter to restrict the amplitude of the voltage that can be stored on the 
capacitor in the integrating network. 


CONCLUSION 


The analysis shown here is built around the saturated transfer character- 
istics of the control loop. In some of the more complex saturated servo 
systems it is doubtful whether the necessary saturated transfer characteris- 
tics even exist as a single valued entity. However, in these systems where it 
does exist, the analysis shows that the dual mode of response exists in some 
degree at any frequency for which 


[G| cos 9 > 1- | (9) 


where |G| is the magnitude of the gain and © the phase margin of the saturated 
transfer characteristic of the loop. Furthermore, when this condition for 
anomalous behavior exists, the input-output characteristic will exhibit a dual 
response when 


Xi > Xo max Sin 0. 
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Here, X; is the input signal, and Xo max is the maximum output amplitude off 
the system at the frequency in question and 6 is again the phase margin of the 
saturated loop transfer characteristic. For these inputs, the output will be 
either the linear value given by 


reo oat 


Let G 


or the saturated value Xx, ae 


The error will likewise be double valued. It will be either the linear 
value, 1 


Sot PCPreT, 


or something between this linear value and \ Xe] = \(X)el cot 6, where (Xi), 
is the critical input value at which dual response begines, as defined above. 
The precise value of \x,| in this latter case depends on the value of |Xz[. 


The predicted anomalous behavior for the two examples used for demonstra- 
tion purposes is closely borne out by the measured results for these systems 
on an analog computer. Figures 12 and 13 show the input-output hysteresis, 
and Figs. 1) and 15 show the corresponding effects on the error. 
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A NUMERICAL METHOD FOR DETERMINING A SYSTEM IMPULSE RESPONSE FROM 
THE TRANSIENT RESPONSE TO ARBITRARY INPUTS 


Norman J. Zabusky 
Hughes Aircraft Company 
Los Angeles, California 


Summary -- The transfer functions of large physical systems must be known if 
they are to be controlled efficiently. Environmental transient test proce- 
dures are described and their advantages presented. Because impulse and step 
functions cannot be obtained from practical test equipment, the problem of 
deducing the impulse response from arbitrary input and response data arises. 


A successful iteration procedure for determining h(T), the impulse 
response, by using the convolution equation, r(t) = dfn e(t--T) h (fT) aT, is 
described. An actual test transient is analyzed. An analog computer tech- 
nique which mechanizes this iteration procedure is presented. The Appendix 
contains a tape multiplication procedure for performing numerical 
convolutions. 


TESTING FOR DYNAMIC CHARACTERISTICS 


Man always tests what he builds in order to determine if his design pro- 
cedures are adequate and to see if his constructions meet desired specifica- 
tions. In recent years dynamic testing of large physical systems, such as 
aircraft, has become increasingly important. One can design a better auto- 
matic pilot with an accurate knowledge of an aircraft's dynamic 
characteristics. 


The best type of test is one which places the object in its natural envi- 
ronment, whether that be pressure, temperature or frequency. For example, if 
a high frequency filter is built with components whose values have been meas- 
ured at- de or’in the low frequency region, then large errors will arise in the 
over-all transfer function. Stray capacitances and inductances which do not 
show up in low frequency measurements cause these errors. An analogous situa- 
tion is found with aircraft. At present they are tested in wind tunnels, and 
‘the information obtained allows a prediction of the transfer functicns 
describing their actual "in-flight" characteristics. The wind tunnel is a 
simulated environment in that an artificial airstream creates forces and 
moments on a small model poised in various positions. As a final check on the 
"in-flight" dynamic characteristics, it is desirable to use the natural envi- 
ronment of "free-flight tests." 


Befcre testing a system, one should have a rough idea of its characteris- 
tics, such as orders of. magnitude to expect in recorded quantities, the order 
of the differential equations describing the system and what nonlinearities to 
expect. This will allow a more effective and meaningful test, and one to 
which data reduction techniques are most amenable. Care must be exercised 
when testing nonstationary systems. For aircraft this means a constant veloc-~ 
ity and altitude during each test, as well as small perturbations around the 
trim or quiescent condition. 
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TEST PROCEDURES 


Dynamic test procedures are broken down into two essential classes: 
steady-state and transient. The first procedure is familiar to most engineers 
and is one in which the input to the system is.excited sinusoidally, and the 
response is measured at different frequencies. This procedure requires excep- 
tional test and measuring equipment. It is very hard (and sometimes impossi- 
ble) to get good sine waves from the elevators of aircraft. Furthermore, 
steady state tests -take a long time and require an object which is stationary 
with respect to its environment. Because of the length of the test, the sta- 
tionary requirement is sometimes impossible to meet. These difficulties are 
the reasons for transient testing. 


In transient tests one excites the input (elevator motion) with some 
function of time and records the response. This gives one a maximum of infor- 
mation in a minimum of time and recording space. For example, the damping 
factor, & » of an underdamped second order system can easily be determined by 
exciting it with a step function. Figure la gives g vs. R, where R is deter- 
mined from the quotient of "amplitude-deviations" of the step function 
response, as shown in the formulas in Fig. lb. I have used this procedure for 
checking the damping of aircraft simulated on analog computers. It gives an 
accuracy of from 1% to 3%. 
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Fig. la - Ws. decrement ratio for a Fig. lb - The step response of a 
second order system. second order system. 


Most often, idealized inputs as steps and impulses cannot be obtained in 
practice. Therefore, we are forced to use what we can get from our program- 
ming equipment. The result is usually an arbitrary function (Fig. 2) which 
cannot be conveniently described mathematically. One chooses the general | 
shape of the arbitrary input so that it reveals characteristics that one is 
most interested in seeing: to avoid nonlinearities, small signals are used3 
to bring out high frequency effects, rapidly varying input functions are 
used. 


It is most desirable tc excite the system with pulse-type and doublet- 
type inputs. An example of the latter is shown in Fig. 2. These minimize the 
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disturbances to the system from trim condition and contain a maximum of infor- 
mation about the natural behavior of the system. Step-type inputs frequently 
hide the natural behavior in the steady-state value reached. Sometimes it is 
necessary to test the system with several inputs. To cover the frequency 
spectrum of interest the duration of the input function is varied, and to 
determine the linear range the amplitude of the input function is varied. 

This will require adjusting sensitivities and ranges of instruments to each 
input. 


P Gr rt) 
INPUT @(t) [ system | RESPONSE r 


ati 


Fig. 2 - Typical input and response curves. 


The lack of an analytical input poses the problem to be solved, namely: 
deduce the impulse response of a system when the input and response are graph- 
ically recorded functions of finite accuracy. Many sources of inaccuracy are 
present during testing, some of which are shown in Fig. 3. In most practical 
applications one can rarely expect to know the full scale quantity being 
measured to more than three significant figures. Usually 3% of full scale is 
considered good. These errors must be accounted for in order to determine a 
reasonably good approximation to the true impulse response. 
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Fig. 3 - A sequence for recording a physical 
quantity. 
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BRIEF HISTORY AND COMMENTS ON NUMERICAL PROCESSES 


Many investigators (see. References) have developed techniques for deter- 
mining the dynamic characteristics of systems from test data. In many cases 
the techniques are for a specific system or a specific form of test data. 

Most deal with linear systems. However, Shimbrot treats systems with nonline- 
arities whose shape is known a priori. Note that a great deal of the work has 
been done by people in the field of aerodynamics. Milliken*s paper presents a 
thorough summary of the work done in the United States (previous to 1951) on 
determining transfer functions from flight-test data. He includes a very com- 
prehensive bibliography. 


Almost all the techniques listed in the bibliography are numerical. Most 
of them deal with the test data in the time domain, using the differential 
equations of motion of the system as a basis for analysis. In the following 
section another technique for determining the impulse response in the time 
domain will be described and illustrated. It assumes linearity and is based 
on the convolution equation. 


Numerical analysis is applied in almost all the above techniques, whether 
it is done consciously or unconsciously. One must be very careful with these 
applications, otherwise erroneous answers will result. Many problems can be 
formulated and solved analytically. However, the transferal of the analytic 
methods to the domain of numerical methods can sometimes lead to false solu- 
tions or difficulties in obtaining the true solutions.t 


The answer (impulse response) obtained by a particular numerical tech- 
nique is a function of the technique and the data. For the problem posed 
above, there is no unique answer. 


CONVOLUTION EQUATION TECHNIQUE FOR DETERMINING THE IMPULSE RESPONSE 


Approximating the Impulse Response 


This technique makes use of the convolution equation? 


*Zabusky , N. J.: "Numerical Methods for Determinging Impulse Responses from 
Recorded Data." Master's thesis in electrical engineering, M.I.T., 1951. 
This paper has also been published by the Servomechanisms Laboratory of 
M.1L ol e 


ty classic example of this arises when solving differential equations by dif- 
ference techniques. If the basic spacing chosen for the difference approxi - 
mation is too large one can get false answers. For example, when solving the 
differential equation of an undamped second order system (linear oscillator), 
too large a spacing will cause the intrcduction of damping and a decrease in 
the natural frequency of oscillation. See the Houbolt approximation in the 
following book: Levy, S. and Kroll, W.: “Errors Introduced by Finite Space 
and Time Increments in Dynamic Response Calculations," Proceedings of the 
First National Congress of Applied Mechanics, pp. 1-9 (1951). 


| fcardner, M. F., and Barnes, J. L.s "Transients in Linear Systems." 
New York, John Wiley and Sons, Inc., pp. 231 and 262 (1950). 
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r(t) = fe @ (t - T) ht (1) af (1) 
1?) 


where: 


Q(t) is the true (or given) input function, 
r(t) is the true (or given) response function, 
hi(T) is the true impulse response of the system. 


The problem of determining the impulse response from the given input and 
response data is one of solving an integral equation. The convolution equa- 
tion is a Volterra Integral Equation of the First Kind. To accomplish this 
ee the true impulse response hi({) is approximated with the function 
His). 


h(T) = x, hy° (fT) + xo ho (1) + <... xe eth) (2) 


In Eq. (2) it is assumed that the functions h;(T) are known. The x; are 
. the unknown constants which characterize the impulse response and are deter- 
mined to "best" satisfy the convolution equation. The choice of the functions 
h; (T) is motivated by physical considerations. The impulse response of any 
system described by linear differential equations with constant coefficients 
can be a sum of polynomials, exponentials, sinusoids and/or any combination of 
the products of these three functions. Hence the basic approximating function 
is the exponential-polynomial product, ~or 


-1 
hi (1) = entil (x4 + XoT+ x3T° + eevee Xp * De (3) 


In the most general case, the exponent Xn+] Ls a complex number. 


Solution of the Problem 


The technique for determining the unknown coefficients, Xi, is numerical. 
The basic numerical process used is that of "performing a convolution." This 
is the process of evaluating the integral of Eq. (1) when © and h are known. 
A "tape multiplication" scheme for accomplishing this is described in the 
Appendix. 


The technique for determining the impulse response employs an iteration 
process. The exponents are modified from one iteration to the next, and at 
each step the coefficients of the exponential-polynomial products are deter- 
mined. With a knowledge of these coefficients, the exponents are modified and 
the iteration process repeated until a "best" solution is obtained. The num= 
ber Csubscript m in Eq. (2)] of basic approximating functions that should be 
used is equal to the order of the system being analyzed. Frequently, the 
physics of the situation will give the required order; if not, one can guess 
at the order. After a few iterations he will determine whether this guess 
yields an accurate fit or whether an additional basic approximating function 
must be added to h({). These points are illustrated by the following 
examples: 
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Pear r(t) and @e(t) are given as graphically recorded functions. Assume 
that hf) has the forn, 


Xoo x Xoot ; 


This second order impulse response is used when one suspects that the unknown 
system is of first order. The additional order provides the flexibility 
required to iterate. In order to construct the h tapes for convolution, x99 
must be known. Usually some a priori information on the system's behavior will 
furnish a value of x99. If no information is available, a procedure which 
decomposes” the given transients and yields a crude step-response is recom- 
mended. Otherwise, one must make a wild guess as to the value of X59 from the 
shape of the given transients. The zero in each subscript of x represents the 
first approximation for each coefficient. The next step is.to solve for x19 
and X39. Substitute Eq. Ory into far (1); 


: Xoo! ts Xone 
ro(t;) = mo J a(t; a e os at + x30) . a(t; = 1 Te oe at. (5) 


The subscript c on r,(t,) designates the "computed" response. This is to dis- 
tinguish it from the "given" response r(t.). The upper limit o can be 
replaced by t; if e(t) = 0 for t<0O. If the need arises, this can be accon- 
plished by a translation of coordinates. 


If we let 
a1 = ils a(t, -7T) a at (6) - 
% : 
t. oh 
853 = ie a(t, = T) ree at (7) 
then 
o(tz) = X104412 + *30%43- (8) 


The "convolution coefficients" a,j and a;, are evaluated numerically as 
described in the Appendix. .If we let tes) = r(t;), we can obtain as many 
equations as there are values of r(ty) given. For example, 


r(ty) = x19 411 + X30 413 
r(to) = X19 421 + X30 423 


r(t3) = X10 437 + X30, 433 (9) 


r(ty) = Xo 4p + *30 2nj 


*Zabusky, op. cit., pp. 127-136. 


In matrix form this is 


ETO 
t =| a.. = oe 
r( ) a4 X50 is X30 (10) 


Now x19 and x39 are to be determined so that ro(tj) will agree most 
closely with the given r(t). Two procedures are available. In the first we 
solve for xj9 and X39 using all the equations in (9) in a least square manner. 
This procedure minimizes the mean-square error “MS, 


Ex “J [r(t) - ro(t) |? at. (11) 
(e) 


Equation (10) is substituted into Eq. (11) and the integration from O->T is 
replaced by a finite summation. If the equations are manipulated, one obtains 
two linear simultaneous equations. 


as | P| et) = |aya| > | ass | . (12) 


Note: j = 1 and 33 [ass | is ani x j matrix whose elements are the convolu- 
tion coefficients and [D| is a diagonal matrix whose elements are 
area-weighting factors. (These correspond to the integration from 
O -T, For Simpson's Rule these diagonal elements are 
Mees i tcceee Uy Le) 


A much easier procedure fits r(t) at as many points as there are unknown 
x's (two in this case). Hence re(t) and r(t) will agree exactly at two points. 


I chose the second procedure because the first involves at least ten times — 
as much numerical work for a fair sampling, and because erroneous solutions may 
arise when solving the simultaneous equations of (12). These false solutions 
result from the accumuiation of truncation and roundoff errors in the addi- 
tional integration and the matrix multiplications. The chief disadvantage of 
the second procedure lies in the fact that xjo9 and x39 are not only functions 
of the assumed x99, but also of the values of t at which the curves were fit- 
ted. The first procedure is worthy of more investigation if accurate high- 
speed computing equipment is available. 


When plotted, rg(t) will most likely differ from r(t), except at the 
points at which they were fitted. If the physical system is truly of first 
crdcer and if the effects of truncation and roundoff errors are negligible, the 
difference is due to a wrong value of x9. To be able to iterate, a new value 
of x9 = X59) must be cbtained. 


x + 
Xo, = Xo + 2 (x) (13) 
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where k is a convergence factor. In the examples chosen a value of k = 1.3 was 
used. When starting the iteration procedure, values of k such that 0.80<k < 
1.30 are acceptable. After a few iterations are performed, k should always be 
chosen less than 0.80. The choice of k will determine whether the iteration 
procedure is stable. With this new value of x9 = Xo,, the above process can be 
repeated. The iteration is continued until x3>Q. If x3 does not appreach 
zero, one must consider four possibilities: 


1. The system is of second or higher order. 
2. The system has a double order pole at Xo (a rarity). 


3. Inaccuracies in the original data and/or numerical truncation and round- 
off errors have introduced extraneous solutions. 


4. The points at which r,(t) was fitted to r(t) have introduced extraneous 
solutions. 


In case 2, r,(t) will agree very closely with r(t) and the problem can be con- 
sidered solved. In cases 1, 3 and }, r,(t) will differ appreciably from r(t); 
and this leads one to use an impulse response of a fourth-order system. 


h(T) = x19 mac (IST S30 X59 eho’ a+ Z708) (14) 
eval) x50 


Actually, a second order system is suspected. The extra two orders in the 
impulse response gives one the flexibility for iteration, i.e., changing the 
exponents. After a number of iterations, one will obtain a better agreement 
between r,(t) and r(t) than was obtained with the second order impulse 
response, Eq. (4). This process of introducing additional orders can be con- 
tinued until the desired agreement is obtained between r,(t) and r(t). 
Unfortunately, no straightforward methods are available for determining when 
extraneous solutions have been introduced. That is why one should have ana 
priori acquaintance with the system and its order. With this knowledge one 


ava Xo 
“This relation is obtained by letting xo = x99 t 5 in the expression for e an 


2 
Se Die series leat eT: (Ge teas 


and all remaining 


If § is very small, then, to a rough approximation, 


(oT) 2 
ie 
terms are negligible. By comparing Eq. (4) (whose first term is normalized) 

x 
30 


with the above equation, one obtains b= cea The convergence factor, k, 


takes into consideration the neglected terms in the series expansion. 
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can start with a good approximation of the impulse response and keep the order 
the same during the iteration. If erroneous solutions are suspected, one 
should decrease the basic spacing used in performing a convolution and repeat 
the procedure. 


To aid the convergence of the iteration process, one can triple the order 
of the approximating impulse response. For example, if one suspects that the 
given data describes a first order system, he can use Eq. (15) rather than Eq. 
(4) to approximate h, (1). 


x ss 
ACI Ne so 6 Oa ee er (15) 
=20 A110 
The values of AS “uO when found, can be used to give a more accu- 


x x 
10 10 . 
rate value of 6 than that obtained from Eq. (13). This procedure increases 
the work involved per iteration but will decrease the number of iterations 
required to obtain a good solution. 


_ The procedures described above can be adapted to systems whose poles are 
complex-conjugate. An iteration procedure based on the minimization of a 
defined error quantity, as fuss will probably prove most advantageous. 
Example Problem 


Figure eat the input and response of an actual physical system. 
From a priori information we write h( oe 


ee eo a, x jen ee + =r (16) 
10 : So 
where 
X9Q = -1.40 
X6Q ® 735-6 (17) . 


The input was convolved numerically with Eq. (16), and the result fitted at 
the values of t aligned with the black arrows (trial 0). Although there are 
only four unknowns in Eq. (16), five equations, and hence five unknowns, are 
indicated by the five black arrows. The extra unknown is a bias which is 
added to r,(t). This takes into consideration unknown biases which arise when 
reading the recorded data. Hence, the general equation used in fitting will 
be of the form 


r(ti) = X19 ai1 + X30 434+ X50 a5 + X70 a7 + XB0- (18) 
where Xg9 is the bias. The trial 0 impulse response takes the form 
hg(T) = 2.66 en1-4OT (1 - 17641) + 4.33 e7 "2357 (1 - 2.737). (19) 
The bias xg, = + l.3. 
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Fig. 4 - The input and response functions obtained from oscillo- 


graph film. 


From this, 6 is computed as (k = 1.30) 


§ 21 = 1.3 (-.76) = -.990 
561 = 1.3 (-2.73) = -3.55 (20) 
Hence Xo, = - 2-39 and x) = - 10.9. The impulse response for the trial x: 
iteration becomes 
= *3]. = cir fib 
ny (P) = xy707239F (a + ma + x51¢ 10.91 (1 + ma (21) 


The five white arrows indicate the points of fitting for trial 1. The follow- 
ing was determined as h,(1): 


h(t) = 1.48 en 2939T (1 + 3587) + 4.76 en2O-9T (1+ 6.uOT). (22) 


The bias for trial 1 xg, = +6.05. Figure 5 shows the error between r(t) and 
ro(t) as a percentage of full scale, that is 


E(t) = ee r(t) - rg(t) + xe) (23) 


Note the appreciable reduction in the error obtained by just one iteration. 
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Fig. 5 = Per cent error as a function of t. 


ANALOG COMPUTER TECHNIQUE FOR DETERMINING THE IMPULSE RESPONSE 


The above procedures are all numerical and as such are adaptable to solu- 
tion by means of desk calculators or digital computers. I would like to 
suggest one of possibly many analog computer techniques that can be used to 
solve this problem. ° 


Figure 6 shows how this technique is mechanized. The given input 6(t) 
and the given response r(t) are produced by analog function generators. The 
input is fed to the analog computer and the response of the computer, r_(t) is 
mixed with the given response r(t) and a proper bias. The result is fed to 
the circuits shown. The output of the top integrator will then be the mean- 
square error, MS. 


T 
Eus = ef Es mere tat (osas)|? dt (2h) 


Oo 
pee the output of the lower integrator will be the mean-absolute-value error, 
MA 


T 
Ea = | ke ~ ro(t) + (sas)| at. (25) 
fe} 


The correct impulse response is obtained when both these quantities are 
absolute minima. For example, if x,, x and the bias are all zero, then the 
error, €, will be only a function of xX, and Xo. This can be depicted by the 
surface shown in Fig. 7. By suitable adjusting potentiometers xX, and x5, the 
point (x1 )min and (Xo)miy will be reached. These two quantities then charac- 
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Fig. 6 - Determination of the impulse response using an analog computer. 


terize the correct impulse response. The logical sequence for varying these 
potentiometers has been described mathematically.” In essence, one seeks the 
point in the error space at which the first. partial derivatives of the error 
are all zero, that is 


a zo , rb ia Bar ct (26) 


This is a necessary condition. 
CONCLUSION 


A numerical iteration procedure for determining the impulse response of 
an overdamped system has been demonstrated. The error of a successful itera- 
tion procedure must converge to a volume in error space which satisfies given 
requirements, and the rate of this convergence must be rapid. Since conver- 
gence depends on the "general shape" and accuracy of the data, no positive 
statements can be made. Future work on this problem should be directed along 
these lines. ; 


It is to be noted that the original data was used in all the above calcu- 
lations. To treat "noisy" data conveniently one must calculate the input | 
auto-correlation function, f 34 (t), and the input-response cross-correlation 
function, 5,(t) from the given data. These are related by the same integral 
equation: 


o0) 
(t) =| h, (t) Pte = od t 
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Fig. 7 - The error in three dimensional space. 
APPENDIX: PERFORMING A CONVOLUTION 


To visualize the process of convolution, first examine the nature of the 
function, @(t -T), designated as the shifting input. The shifting nature can 
be visualized by considering the simple example in Fig. 8 (A, Band C). Here 
eT) = +f is a straight line of slope +1. 0@( - 7) is obtained by replacing 
+ { by-T, and the result is a straight line of slope-1. @(-T) is a 
mirror image of 6(f) around {\= 0. If we replace ( -T) by (t -T), then 
e(t -T) = t -T is a straight line of slope - 1, shifted in the positive 
direction by an amountf =t. A completely analogous discussion holds for any 
arbitrary function as shown in Fig. 8 (D, E and F). As t takes on increasing 


values, the mirror image (or folded) curve shifts along the abcissa in the 
positive T direction. 


To perform a convolution, one shifts the input function by an amount t, 
multiplies its ordinates with corresponding ordinates of the impulse response 
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O(t-2)= t-7 


Fig. 8 - The shifting input. 


and obtains a product curve, P(t -T), as shown in Fig. 9. The area under 
the product curve is the required convolution. 


Numerical convolution is accomplished easily with convolution tapes. 
They systematize the calculations for desk calculating machines and minimize 
the writing and thought required. The "h" tape is constructed by recording in 
boxes the ordinates of the h curve at an equal spacing 4. (See Fig. 10.) 
This tape starts at T = 0 and runs toward the right, just as the impulse 
- response starts at T = 0 and runs toward the right. The 6 tape starts at zero 
and runs toward the left. In each box of this tape one records an ordinate of 
Q(T) each weighted ur 2 a factor depending upon the mumerical integration scheme 
being used. Simpson Rule is Sean in the @ tape of Fig. 10. 


To evaluate r,(8), one shifts the @ tape to the right by 8& units as 
shown in Fig. 10. Then one multiplies the values of the ordinate h by the 
weighted value of @ in the box immediately below, and adds the product to the 
previous products. The quantities in the second row of the @ tape are used 
only for the last multiplication in any convolution. Hence 


re(@d) #8 [0(o) h (BS) +u e(S) h (TS) +2 (2B) h (6S) +h OC) 
th (St) + 2 O(LN) bh (UES) + 4 0(5ts) h (30N)-+ 2 0(6D) h (205) +h O(7E) 
h (%) + 0(88) h (0). (27) 


The operations of multiplication and accumulation are accomplished conveni- 
ently with a desk calculator. . 


“7 abusky, Ope cite, PPpe 78-87. 
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IMPULSE RESPONSE = h(7) 
Tv 


t—+ 


PRODUCT CURVE = h(%)6(t-t) = P(t-2) 


Fig. 9 = Visualizing a convolution. 


The size of the interval & is determined by the shape of the product 
curve, P (t -{*). &S should be small enough so that a parabolic approximation 
to any segment of the product curve of 2 ( units in length will be an accu- 
rate approximation. 
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Fig. 10 = Convolution by tape multiplication. 
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SYNTHESIS OF FEEDBACK CONTROL SYSTEMS WITH A MINIMUM LEAD 
FOR A SPECIFIED PERFORMANCE” 


George S. Axelby 
Westinghouse’ Electric Corporation 
Baltimore, Maryland 


Summary -- Compensation is usually needed in a high performance feedback 
control system to realize a desired frequency response. Direct synthesis pro- 
cedures to determine the necessary compensation are described. In particular, 
consideration is given to "minimum lead" systems which have an optimum gain. 


SS 
_The desired performance of a feedback control system may be specified 
approximately in different ways. In some applications, the transient response 


of a system to a particular input may be prescribed, and in other systems the 
frequency response may be specified. 


Of course, these criteria are related to each other in a given system, 
but system performance may be conveniently specified approximately in terms of 
the frequency response, as summarized in Fig. 1: The open loop and closed 
loop frequency responses and transient responses of the system are indicated, 
as well as the approximate relations between them which hold for a large class 
of systems. . 


Using these relationships, the band-pass, G)., the position constant, 
velocity constant, the minimum phase margin or the phase margin, the frequency 
of minimum phase margin and the gain at particular frequencies may be deter- 
mined approximately from desired performance specifications. Of course, the 
requirements may conflict if all variables are specified. In this case, the 
usual design compromises are necessary. 


In the following discussions, it will be assumed that the desired system 
performance can be translated into the frequency response criteria shown in 
Fig- 1, and synthesis techniques will be established in terms of the frequency 
response. 


METHODS OF SYNTHESIS 


Most of the synthesis methods currently in use employ cut and try methods 
to realize performance specifications. This is not very efficient and fre- 
quently the most feasible arrangement of parameters is not found because of 
the time absorbing details. 


Often it is difficult to obtain an over-all perspective of system perfor- 
mance. ‘Special computing aids and charts are usually used to speed synthesis, 


“This paper was presented at the 3rd Southwest Conference on Feedback Control 
Systems in Dallas, Texas, April 28, 1955. The conference was sponsored by 
the Department of Electrical Engineering, Southern Methodist University in 
cooperation with the North Texas Section, A.1.E.E., the Professional Group on 
Automatic Control of the Dallas-Ft. Worth Section, IRE. 
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OPEN LOOP. FREQUENCY RESPONSE 


Ke lopprox.) 
ah lege \ 


Lales|*|-< 


Gyr:G,'G,°G, 


4 K\(s/wotl) 5 b ft 
T (s/wstl) sls/w, +1) 


APPROXIMATE PERFORMANCE 
CRITERJA 
(A REASONABLY STABLE SYSTEM 


1S ASSUMED,) 
=2 


FOR AN INPUT, F,, WITH A DOMINANT 
FREQUENCY COMPONENT, Wy, OF 
MAGNITUDE, |R|, THE ERROR 
MAGNITUDE 


lel S IRles, 
ISlw, 


FOR AN INPUT, r,, CHANGING 
SMOOTHLY, THE TIME ERROR, 


TRANSIENT RESPONSE 
1 CLOSED LOOP 


1.0 ae 1.0 


M 63 
\c/R| Some 
OUT PUT,C 
(UNIT STEP INPUT) 
0 Wp 
J 
-90° ERROR,E 
(UNIT RAMP INPUT) 
Me 
-180° 


Fig. 1 - Approximate performance criteria from the frequency response. 


but they are not always available, and they do not always show basic relation- 
ships between system parameters which could improve design. However, charts 
an2 computing aids have an important part in system design, particularly in 
analysis and in checking the quciker, more direct, but less accurate method to 
ve described. Although it is sometimes inaccurate (there is a warning when 
tnat inaccuracy is present, however), it is a powerful addition to the many 
analysis and synthesis techniques now in use. It can be used to obtain gen- 
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éral relationships between system parameters which uni.quely optimize the 
system to meet design specifications. 


From simple, approximate equations, the necessary system compensation 
and transfer function may be directly determined, once the pertormance is spe- 
cified, in terms of system gain, band-pass and degree of stability. This 
synthesis technique will be referred to as "asymptotic synthesis." 


"MINIMUM LEAD" SYSTEMS 


Because an infinite variety of feedback control systems could be consid- 
ered, the discussion will be limited to a type of system, described below, 
common to many high performance servomechanisms. 


It may be noted that high performance systems often have large integral 
‘networks or double integrators, a high gain at low frequencies, and a minimum 


phase angle, ©, which.occurs near the crossover frequency.. This is illus- 
trated in Fig. 2. 


LOG MAGNITUDE OF 
OPEN LOOP TRANSFER 
FUNCTION, G 


PHASE OF G 


Fig. 2 - Asymptotic gain and phase angle of G, 


The minimum phase angle occurs at a frequency @,,, below the open loop 
crossover frequency, 7,3; Gz, the open loop gain at @m, equals ©,/#,,, assuming 
that slope of the open loop magnitude is -1 where the open loop magnitude is 
unity. For a given M,, the maximum magnitude of the closed loop response, it 
would be desirable to have the frequency response locus of G, when plotted on 
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a Nichols chart, tangent to the given M circle at its maximum excursion from 
180 - This is indicated in Fig. 3. The gain at this point is an optimum for 
a given Mp because a higher or lower gain will produce a larger Mp. 


In the design, this is accomplished by placing wa at the desired tangént 
point. Thus, for a given Mj, @, should have the following relationship to >: 


OD, =&%, cos Dm = 


wo 
i} 


LOOP TRANSFER 
LOCUS OF G 


DESIRED Mp Mp sobs 
LOCUS OF MAXIMUM PHASE MARGINS 


FOR A CONSTANT Mp 
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1) 


We 


ig 2g 
MMM EL cosdm Ym 


Mp? Tee 


Fig. 3 - Nichols chart of minimum lead system. 


Such an arrangement provides the desired stability with a minimum amount 
of effective lead (as denoted by the -1, or -6 db/octave, slope through 4.) 
for a particular crossover frequency. This often economizes the amount of 
compensation needed for stability; it helps reduce noise within the system and 
optimizes the system gain. Furthermore, the transient settling time of such a 
system will be close to a minimum for a given band-pass. 


A system with this design will be referred to as a "minimum lead" system. 
ASYMPTOTIC SYNTHESIS METHODS 


The method of asymptotic synthesis is relatively simple. The desired 
system transfer function, which appears to meet the most severe of the many 
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possible specifications, is sketched roughly to obtain the asymptotic form of 
the transfer function. Usually several of the corner frequencies will not be 
known, and it is desired that these corner frequencies be determined to pro- 
vide the necessary, if not optimum, gain band-pass and degree of stability. 


Three simple basic equations, and possibly others, are written. The 
asymptotic sketch of log gain vs. log frequency (log-log coordinates) is a 
valuable guide when writing the equations and when making the final check of 
the calculated system. These algebraic equations are written to relate: 


1. the minimum phase angle of the open loop near the crossover frequency 
to the known and unknown corner frequencies of the system, 


2. the system corner frequencies in a manner which will establish a "min- 
imum lead" system, 


3. the desired open loop gain to the known and unknown corner frequencies 
and to the system crossover frequency. 


THE BASIC EQUATIONS 
These sets of equations contain common unknowns, and the equations may 


be solved simultaneously. The method of deriving these equations will now be 
discussed, and applications of them will be given. 


Equations relating the phase angle of the open loop near the: crossover 
frequency (or between any corner frequencies) may be written as follows: 


Consider the following open loop transfer function: 


Geter s)s (sr rss) 
~ (1+ 73s) (1 + Ts) 


z = Sareea alee 


J@, : (1+ = (1+ poe ee 


where: 4), ete healed at which phase is desired, and 


As wae 3, a) = corner frequencies. 


From G the phase shift, P,, may be written: 
je : 
c 


P. = are tan GE) + arc tan G+ «se ~ are tan @) 


—are tan (=)... eee 


Each term on the right may be represented by an infinite series. 
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It would seem that the use of an infinite series in place of the arc tan 
of a number would infinitely complicate the problem of finding the phase 
shift, FP ., of G s - A very good approximation .can be obtained, however, by 

J~C 

using the first term of the series when#,/#, £ 1, and. by ee the first two 
terms of the series when4%,/a) > 1. Thus; are tan (flr) @),/t),, when 
@),/J; <— 1 (corner frequency above ©.)3 = W/2 -G/e},, whend ae; > 1 (cor- 
ner aa below %). (If a de ae nonianto quadratic should exist, the 
approximate phase shift at @, would be: -2§ G4/@,, Wy > «4; or 

-1T + 2g 44, CO), »®, <4, >) nots that the frequency ratio (of corner frequency 
and Seek es BS which phase is desired) is always less than 1 in ‘the 
agrreximation. 


This approximation is especially good if the frequency ratio in the 
approximation is equal to or less than 0.5. In other words, the approximation 
‘ hoids particularly well when the frequency at which the phase shift is desired 
is at least 2:1 from the nearest corner frequency. The crossover frequency is 
usually in such a position, and it is in this region that the phase shift is 
most important. Frequently the phase shift is specified in this region 
because it is here that the phase shift determines the degree of system 
stability. | : 


The error due to the arc tan approximation is shown in Fig. 4. Correc- 
tion factors could be used to minimize the error, but errors tend to cancel, 
and the method is usually accurate enough without them. Correction factors 
would only add complexity, and the advantages of this synthesis method would 
be decreased. 
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Fig. 4 - Comparison between actual and approximate phase shift for a single 
asymptotic corner as computed by conventional and approximate 


methods. 
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Thus, an equation for the phase shift can be reduced from a summation of 
trigonometric terms to algebraic terms. As shown in Fig. 5, an expression for 
the phase angle, Bc» at We can be written from inspection of the asymptotic 
diagram, and 


figs (le Ol) + eL, = pe 2 &) 
Bo 2h Geter EG, 2 7 W ~ G3 i 
Note that: 
1. - 1/2 occurs because of the -1 slope, 


2. a term appears for each corner frequency, 


3. each term is multiplied by plus or minus the Change of slope at the cor- 
ner frequency, pee 


4. the corners.may or may not be known. 


Thus, B. can be expressed as an algebraic equation 


eee eee ao Q 
Po 2 +E - Uy ~ Bs 2D, 


This is the first of the three basic system equations. 


The second basic equation is that which relates the system corner fre- 
quencies in a manner which will establish a "minimum lead" system. Of course, 
this may not always be required; other specifications may take precedence. 
This specification may unduly complicate the problem, or it may produce redun- 
dant equations if too many variables are specified. 


This equation is found by differentiating the first equation for B 


(an 
Thus, 


Although Yo, Q;, or ©, may be determined from the last equations in terms of 
other frequencies, it is usually more convenient, for later substitution, to 
separate the corner frequencies above and below (J), as shown. 


_ The third basic equation relates the desired open loop: gains to corner 
frequencies, known and unknown. 


This can be developed best by referring to Fig. 6. 


Consider the gain between the frequencies any and @5, On log-log coordi- 
nates the gain change is a straight line with a Slope of “231. 


62 


FORM FOR FORM FOR CORNER 
CORNER FREQUENCIES<We FREQUENCIES >we 


In GAIN 


BY INSPECTION OF THE ASYMPTOTIC PLOT, THE PHASE ANGLE, Be, AT 
We IS APPROXIMATELY: 


Be. |S IN RADIANS 


NOTE THAT CHANGE OF SLOPE AT CORNER FREQUENCY. 
PRECEDES THE ARC TAN TERM. 


PHASE MARGIN, y, at we = 180° + 22 A, 


NOTE: FREQUENCY SCALE IS IN LN. RADIANS TO RELATE CORNER 
FREQUENCIES TO TIME CONSTANTS, AND THAT —l, —-2, SLOPES &-6, 
—|2, db/OCTAVE 


Fig. 5 - Approximate phase calculations. 
Thus, by the definition of a -2 slope, 
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Fig. 6 = Relationships between gain and phase. 


For ¥ n slope between 4, and ®o, 


Thus, by inspection of the asymptotic diagram in Fig.°6, the velocity 
constant of the system, 
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This algebraic equation must always hold true for the configuration shown 
in Fig. 6. It also serves as a check of the asymptotic drawing. ; 


In a similar manner K, We oe 
s er spore Vana 


K,» 4 &, 


and similarly, the asymptotic gain may be specified at any frequency as a 
function of the corner frequencies. 


To illustrate the use of these equations » consider a typical system in 
the following problem: 


Suppose that a sinusoidal input with a frequency of 1 radian/second and 
an amplitude of 1° peak is applied to a position system. It is desired that 
the peak sinusoidal error be less than 0.05° peak. It is also required that 
the lowest possible system gain and band-pass be used and that the frequency 
response peak magnitude be equal to or less than 1.3. It is known that a sin- 
gle integration exists, and that a corner frequency occurs at 180 rad/sec. It 
is desired to attenuate frequencies above the crossover as soon as possible. 


The form of the open loop function that will meet the requirements is 
shown in Figs. 2 and 3 except that @, will eqaul 180 rad/sec, and it will'con- 
sist of a single corner frequency. 1 will be placed at @/= 1 rad/sec, and 
the asymptotic gain at 0] = 1 (at the corner frequency) will be 
V2 (1)/.05) = 397 = 00. This is also the velocity constant. 


From the specifications, it will be necessary to have a Minimum lead" 
system as previously defined. . 


This means that the minimum phase shift in the, system must occur below 4 
at Ae where 


4), = W, cos f, 
G), is the frequency at which the minimum phase B.. occurs, and 
Dm = 180° + &. 
The basic calculations, the three approximate equations, are derived in 
Fig. 7. They are also shown below: 


1. Phase shift ata: 


2. Condition that makesP,, a minimum: 
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2. Relation between Ky, and corner frequencies: 
A 
kK, = a), me 
mak 


Three other basic equations, derived in the Appendix, complete the mathemati- 


cal description of the system. For a minimum lead system, these relationships 
are exact. . 


sin Bra = 1/My 
Gh =, cos By 


z, ul I | ! 


2 am (2-H) -8m(aS + igo) 


(THE ANGLE AT Wp, NOT we, IS SPECIFIED) 


(NOTE THAT THE FREQUENCIES ABOVE AND BELOW 
Wc ARE SEPARATED) 


THE THIRD EQUATION: 


(st) 


THESE EQUATIONS ARE GENERAL, AND WITH STUDY, USEFUL 
INFORMATION CAN BE OBTAINED FROM THEM. THERE ARE 

9 UNKNOWNS; 6 EQUATIONS. IF 3 PARAMETERS ARE SPECIFIED 
THE REMAINING 6 WILL BE DETERMINED. 


Fig. 7 - Asymptotic synthesis equations. . 
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These equations are quite general and with study useful design formulas can be 
obtained from them. There are nine unknowns, six equations for the system 
shown in Fig. 2. If any three parameters are specified, the remaining six can 
be determined by solving the equaticns simultaneously. This may be done 
rapidly and the results are surprisingly accurate if 


2 £0, Zz O3/2, 


For the problem in question, three variables have been determined from 
the specifications: 


@ =1 My = 1.3 Ky, = 00 
These values have been substituted into the three equations and the 
simultaneous solution is given in Fig. 8. From the solution, 


a 9.9 rad/sec 


OE 
@, = 0.5 rad/sec 
Om = 25.9 rad/sec 
WJ. = 129.5 rad/sec 


and the open loop transfer function, G(s) becomes: 


Pee OOG tay 
s(s +1) G35 +1) (a5 +2) 


DISCUSSION 


1 Accuracy of Calculations 


The equations in Fig. 8 constitute the complete calculations. The six 
unknowns were uniquely determined. Note that 24) = 51.8,@,, = 12.95, 
M3 = 129.5 and, = 9.9. This shows that the condition 24 £4, = @4/2 has 
been met, and that the approximate phase equation (1) should have been quite 
satisfactory. 


Only a small error in these first calculations should be expected. Cer- 
tainly the calculations are more accurate than the system parameters can be 
measured. Nevertheless, no approximate method should be trusted. Perhaps no 
one solution should be judged as final until it has been checked by another 
method. This solution was checked for possible revision by graphical methods. 


The corner frequency-gain relationships were checked with an asymptotic 
plot. This is shown in Fig. 9. 


The Mp and @,, were checked by substituting s = jWJ in the transfer func- 

tion and calculating the actual G - This was plotted on a Nichols chart.. 
jo 

The results are reproduced in Fig. 10. This was the first calculation. 
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FOR THE PROBLEM IN QUESTION, 3 PARAMETERS ARE SPECIFIED: 
w=! Mpzl.3  Ky=400 
THEN FROM FIGURE 7, 
$m= 50.3°=0,.88 RADIANS; Am:-129.7°:-2.26 RADIANS ; 


Wm 


(w2-I)(wo) (G00) (a) _ M2we-l) | (We~-w)) (we) 


TaD Ppt np ter ~(0.64)(400) “ (400) (0.64) 


0-688 = (064400) 


(0.688) (0.64) (400 Z 
2 


2 
Wo 


_ it/1¥4 (88) 
= aD ee 


= 9.9 RAD/SEC 


We = 40.5 RAD/SEC | EQ, (I) 
Wm= 25.9 RAD/SEC | Wm=0.64 We. 


FROM Eq.(b), =o,-—-. oF W 3 129.5 RAD/SEC 


AND THE SYSTEM OPEN Loop FUNCTION } g 2 400 


Fig. 8 - Solution of equations in Pape 7s 
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Fig. 9 = Asymptotic plot of synthesized system. 
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180 


WHERE s=jw 


Fig. 10 = Nichols chart of synthesized system. 


2. Limitations in Using the Method 


Unfortunately, all of the uses, advantages and disadvantages of this syn- 
thesis method cannot be included in this brief discussion. Examples have 
served to-illustrate its advantages and its scope. 


Disadvantages of this method exist largely when it is improperly used. 
For instance, a complex system may yield an excessive number of equations, and 
the simultaneous solution may result in an equation of high order. Considera- 
ble time could be spent solving the equation for its roots, only one of which 
would be useful. Possibly a trial and error solution could be made more 
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quickly, but it would have to satisfy very closely the original simultaneous 
equations, written from the asymptotic diagram. Actually, the equations can 
usually be solved rather quickly because the allowable range of each corner 
frequency can be observed from the asymptotic diagrams, and considerable sim- 
plification can be made if some of the variables are known to be much larger 
than others, or if the range of variation is known. Thus, even for complex 
systems, asymptotic synthesis can be used advantageously. 


Another apparent disadvantage in the method is that solutions may yield 
negative frequencies, complex frequencies which were not expected or corner 
frequencies which change the shape of the original asymptotic diagram. These 
solutions would be incorrect, but this is not usually a disadvantage or a 
detriment in the methods; they merely indicate that the proposed system cannot 
meet the design specification, or that more restrictions have been placed on 
the system than the number of variables allow. It is necessary in these cases 
to reshape the asymptotic diagrams, reduce the severity of design specifica- 
tions or compromise the requirements to prevent conflicting or redundant 
equations. 


These modifications can be made quickly using asymptotic synthesis proce- 
dures, an over-all perspective of the possible system performance can be 
retained and the system design can be readily established. 


3. Usefulness of the Method 


When properly used and when its limitations are realized, the method of 
asymptotic synthesis can be a powerful tool to develop rapidly over-all rela- 
tionships between system parameters and to establish a framework for optimun 
system design. It does not replace more detailed or more accurate studies; 
rather it supplements them. It reduces the problem of loop design to rather 
simple algebra which yields fairly accurate solutions directly with a minimum 
of trial and error. 


APPENDIX 


Derivation of the Locus of Maximum 
Phase Margin for a Given M, on a 
Nichols Chart 


N= (pea 
Gc = al € +5? 


Pz-T + Bp 


(Refer to Fig. 3 for symbols) 
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PERFORMANCE OF DRIV MEMBERS IN FEEDBACK CONTROL SYSTEMS” 


F. M. Bailey 
General Electric Company 
Schenectady, New York 


Summary -- Limitations of system performance due to a mechanical structure are 
discussed, and practical methods of improving the performance of fire control 
systems with this limitation are illustrated. 


The history of high performance feedback control systems has been closely 
associated with the development of naval gun power drives. Progress in this 
field had its start with the development of early thyratron, hydraulic and 
then amplidyne systems, and history has shown a preference for each depending 
upon the degree to which the components have been perfected at any given time. 


Closely allied with this work has been the progressive development of 
fire control systems, an important element of which is the director. The 
feedback control problems of the director are in many ways similar to that of 
the gun mount drive, but it is the significant differences which are to be 
discussed in this paper. 


PERFORMANCE REQUIREMENTS 


In discussing performance we must first examine the requirements which a 
good system must meet. In reviewing this, the discussion will be restricted 
to the use of opened and closed loop frequency response diagrams which form 
the basis of system response criteria. 


A simplified closed loop feedback system is shown in Fig. 1. It contains 
an amplifier, a power boost element and a driving element such as an electric 
motor or hydraulic piston. The driving unit actuates the load and aL eriirt . 
addition, connected mechanically to a selsyn control transformer which gener- 
ates a difference signal between load position and order position. This 
signal energizes the servo amplifier. 


Frequency response measurements may be made on the system by cutting the | 
loop at the point (X) and measuring output shaft motion with an oscilloscope 
while introducing a constant amplitude variable frequency signal at the input. 
The open loop diagram thus obtained is shown in Fig. 2. For constant input, 
the output will drop off at the rate of 6 db per octave until the first time 
constant is reached (usually at the motor constant), at which point the curve 
will begin to attenuate at 12 db per octave and so on at an increased rate of 
6 db per octave for each successive time constant. If the loop is closed with 
a gain represented by the magnitude A, the zero db axis will cross the curve 


“This paper was presented at the 3rd Southwest Conference on Feedback Control 
Systems in Dallas, Texas, April 28, 1955. The conference was sponsored by 
the Department of Electrical Engineering, Southern Methodist University in 
cooperation with the North Texas Section, A.I.E.E., the Professional Group on 
Automatic Control of the Dallas-Ft. Worth section, IRE. 
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Fig. 1 - Servo control loop. 
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at 6 db per octave and the system will be stable. The corresponding closed 
loop frequency response is shown in Fig. 2. 


If the gain is increased as represented by magnitude B, the zero db axis 
will cross in a region close to 12 db per octave and the corresponding closed 
loop curve in Fig. 2 will show considerable resonance. If this process is 
extended by increasing the gain to a magnitude C, it is clear that the system 
will become unstable and break into oscillation. This then is basically the 
servo problem -- to raise the gain sufficiently to extend the frequency 
response of the closed loop diagram and at the same time prevent the approach 
of instability in the enclosed feedback loop. 


TYPICAL PERFORMANCE 


With this brief background, Fig. 3 can be examined to obtain a feeling 
for the type of performance which can be obtained with currently available 
servo system components. 


This figure is a composite of the frequency response and phase shift 
curves of three of our major power drive systems. Curve A is an 8 hp unit 
capable of delivering 0 hp peaks and is typical of equipment currently being 
used in the field. It will be noted that the frequency peak rises to nearly 
1.6 and that at 24 cycles, there is approximately 50° of phase shift. Curve C 
is the response for a more recent production unit having a 25 hp rating with 
capabilities of 00 hp peak. This is a cascaded generator system and it will 
be noted that the phase shift at 24 cycles is also 50°. On a test stand, the 
power drive by itself is capable of performing with a phase shift of 30°. 
Curve B shows the performance of an ignitron drive in the same power range as 
system number one. In this case, the phase shift at 23 cycles has been 
brought down to 25°, 


The dotted line connecting the two humps in the response curve represents 
the performance achieved with the unit on a test stand in which a flywheel is 
used to simulate the load (instead of gearing and mount inertia). It will be 
noted that the mount resonance shows up in reality as an antiresonance. This 
is because the position loop is-closed around the drive motor by means of 
additional instrument gearing and the gun mount is coupled through the gearing 
to the servo system. Experience has shown that closing the loop around a set 
of power gearing will lead to stability problems which, seriously hamper the 
performance of the over-all equipment. The performance shown by these curves 
can be achieved only by the device of utilizing additional instrument gearing. 
It is also clear that in searching for increased servo performance, a limita- 
tion, even with this arrangement, is imposed by the mechanical structure 
utilized in the moving assembly. With equipment of reasonable mechanical con- 
figuration, and especially that hampered by the necessary emphasis on minimum 
weight, it seems unlikely that any appreciable system improvement can be 
achieved through servomechanism component improvement alone. 


THE DIRECTOR PROBLEM 
It has been pointed out that in gun drives it is virtually impossible to 
stabilize a high performance unit by taking the response from the one speed 
shaft, and that it is necessary to utilize instrument gearing from the motor 
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shaft for this purpose. A director drive, on the other hand, whose function 
is to force the radar antenna to track a target, must of necessity obtain its 
information from one speed motion since it is at this point that the antenna 
is connected. In this case, the power gearing is included in the closed loop 
and the additional time constants contribute to the stability problem. Reso- 
nant peaks for conventional units of this type come dangerously close to two 
and the phase shift involved may be well over 100° at resonance. Such per- 
formance is tolerable for a chain type of director system; but with continuing 
interest in the super-position of director and firing equipment, the resonant 
problems of the disturbed sight type of configuration become very critical. 
Here, in addition to the complexity of the gearing problem is added the addi- 
tional complexity of an antenna servo all cascaded within the same control 
loop. The elimination of one or the other seems to be necessary in order to 
accomplish the end purpose. 


Presently there are two methods of doing this. The first is to build an 
unwound director system in which a mechanically driven platform is mounted 
above the gun turret and is positioned by unwinding gearing. The second is to 
provide a power drive system in which there is no gearing. A discussion of 
the second system is within the scope of this paper. 


Figure ) shows a schematic diagram of the arrangement utilizing a one 
speed motor for driving either a director or gun platform as compared with a 


DIRECT DRIVE MOTOR 


[i 


PLATFORM 
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ARMATURE 


GEARED MOTOR MOTOR OPERATES AT LOAD 
SPEED, NO GEARING 


Fig. 4 
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Fig. 5 


geared train system for doing the same thing. Figure 5 shows a photograph of 
an electric motor built for this purpose weighing 700 lbs. and providing an 
output torque of 1000 ft. lbs. 


The frequency response of this unit is shown in Fig. 6. As far as can be 
ascertained, utilization of such a piece of equipment can make a disturbed 
line of sight director system successful. Such a unit employing a conven- 
tional style motor has been set up at the Schenectady County Airport and the 
tracking and computing operation accomplished utilizing a gyro lead angle 
computer. 


It is realized that the weight considerations are entirely unsatisfactory 
for airborne equipment, and for this reason we have initiated a program to 
accomplish the same functional operation with hydraulic equipment. 


Figure 7 presents a picture of a five piston hydraulic actuator which can 
be used as a one speed element in a power drive system. The building of such 
an equipment has been very recent and performance data is not available 
although the operation of similar hydraulic equipment leads us to be very 
optimistic. 


79 


9 °STd 
QNOOSS Yad SNVIOVY 


| 
a a nx © n © ry - 


AE Hee 


et a, PETE 
Eel ee 2 Ee SEE 
eee | eee 


2 


re 


ao 
PTE 
Rept 4 


i Hall 


E 
| 


hi 
Hh 
f 


tl 
i HH 


ae 


aia 


ca 


i +4 14. 
alee fH SES EE 
¥ F Th 1 
rhs i a Cy 


SSS rere SS 
ical Suid Git etn ith SS eae 
i 


80 


NEW COMPONENTS FOR HIGH PERFORMANCE SYSTEMS 


é The previous discussion has demonstrated the benefits obtained from 
introducing simplicity in the feedback systems. The contributions in this 
area can be best provided by making components which, when combined, can form 
the simplest possible servo configuration. Since performance can be increased 


Fig. 7 - Five piston hydraulic rotary actuator. 


by including a rigidly mounted tachometer, it is essential that we build actu- 
ating equipment such as servo motors with integrally mounted tachometers. For 
practical purposes nearly all servos require a brake, and the construction of 
a combination motor, tachometer and brake package is of advantage. 


In Fig. 8 are shown two electric packages -- one a 2 hp, and the other a 
4 hp unit capable of delivering 25 ft. lbs. and 10.ft. lbs. of torque, respec= 


tively. A second step is to provide a power amplification unit which has a 


minimum number of time constants to be used to actuate the motor. In these 
horsepower ranges, such a system seems to be a generator since it offers the 
servo designer one time constant instead of two found in a normal amplidyne. 
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Fig. 8 


Fig. 9 


For higher power ranges, two stages of generation may be necessary in order to 
get adequate power gain, but it has been demonstrated that this arrangement is 
preferable as compared to compensating for the time constants involved in 
large size amplidynes. A sample of such a generator is shown in Bae Fe A 

A third approach is to provide a series of one speed motors which can be 
incorporated into servos to eliminate the gear train. Work covered to date in 
this area has been discussed in the previous section. 


There is one other area in which hope for improved performance lies and 
that is in the operation of the position servo as a combination velocity and 
position system. This arrangement supplies a gun order signal which includes 
both a position signal and a velocity signal. The fundamental scheme to be 
considered centers around the fact that a velocity system has 90° more permis- 
sible phase shift within the closed loop and hence can tolerate more time 
constants than a position system. Therefore, if our feedback system is built 
up on the basis of a velocity system and sufficient gain is introduced to 
maintain dynamic shaft correspondence, there will be little difference y 
observed in the performance of the system except for drift as compared to a 
position system. The drift problem can then be corrected by means of a low 
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gain stable position tie. Such equipment has been built and demonstrated in 
the case of naval directors and experimental work on army tank stabilization 
units. To design a servo of this type, there seems to be a need for an inte- 
gral synchro/tachometer combination to simplify the order signal transmission. 
Problems of linearity and voltage gradient must still be examined before ade- 
quate equipment of this nature can be built. 


CONCLUSION 


A review of the type of performance which can be expected in modern feed- 
back systems indicates that progress will be made by improving components. In 
this paper much emphasis has been placed upon the activities in the field of 
naval equipment where horsepower demands are large and where size and weight 
are not as dominating as in the case of aircraft equipment. We are not over- 
looking the aircraft problem, however, and in certain areas have already made 
attempts to develop the hydraulic analogs which will bring lightness of weight 
to the same applications which are apparent in naval equipment. The develop- 
ment of improved hydraulic components holds much promise that we will be able 
to do this. 
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AUTOMATION AS THE ENGINEER SEES IT” 


W. R. G. Baker 
General Electric Company 
Syracuse, New York 


Summary -- Problems of automation are discussed with particular emphasis on 
their importance to society and to the engineer. 


What makes automation such a fascinating topic to engineer or social sci- 
entist, manager or labor leader, is that each views automation from the 
perspective of his own experience. This perspective, sometimes creates dis- 
tortion. However you think of automation depends entirely upon your point of 
view; to my mind, the answer lies considerably to the right of center. 


How much my outlook toward automation is colored by my engineering back- 
ground is a matter for conjecture. I sometimes find myself agreeing with the 
economist rather than with the technical man. In a recent Harper's article, 
economist Peter Drucker spells out a number of truths about automation. | 
First, he states that automation is something that is here now, not sométhing 
that is approaching; secondly, that automation is not just the pushbutton fac- 
tory. Mr. Drucker then states that there are three basic principles that make 
up the logic of automation: economic activity as a process}:-the pattern, 
order or form behind changing economic phenomena; and finally, self-regulating 
control. He closes his article with the statement, "Automation is not gadg- 
eteering, it is not even engineering; it is a concept of the structure and 
order of economic life, the design of its basic patterns integrated into a 
harmonious, balanced and organic.whole.” 


In technical terms, and limiting the definition simply to the mechanical 
aspects, automation is the addition of the tools of programming and computing 
to the well used tools of mechanization -- mechanics, hydraulics and pneumat- 
ics. If this be automation, how then can it also be "a concept of the 
structure and order of economic life"? I think the apparent disparity in 
viewpoints can be reconciled. To the engineer, the machine is predictable, 
because of the physical laws that govern its actions. To the engineer, soci- 
ety is often unpredictable because action and reaction are sometimes not t 

immediately apparent or may be obscured by the reaction to a combination of 
forces. 


To the economist, society is not unpredictable if all of the economic 
forces at work can be discovered and identified. And if the self-regulating 
factors influencing society are recognized, then "feedback" becomes one of the 
essential elements in the structure... 


To the engineer, then, automation relates merely to the machine; but the 
engineer can see that the utilization of the machine may have far reaching 
economic effects. This may become a disturbing factor in the engineer's 


“This paper was presented at the Centennial Symposium, School of Engineering, 
Michigan State College, East Lansing, Michigan, on May 13, 1955. 
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thinking, for he is not certain that what he is building is beneficial for or 
detrimental to society. It becomes advisable for the engineer to broaden his 
thinking, so that he may see the "animal" as a whole, and not as a mere 
machine. Automatic machinery, or the automatic factory or office, when seen 
in proper prespective, becomes one result of automation, in this broader con- 
cept, and not a primary cause of automation. 


Automation is not a thing of the future; it is here avid is rapidly 
increasing in scope. We must talk of it in terms of the present. 


There are instances where consumer products plants are now turning out 
appliances at the rate of one every 30 seconds, with 2500 machines strung over 
25 miles of conveyors. In data processing, large computers have taken over 
the work af handling payrolls for thousands of employees. As many as three 
million accounts of a national insurance company are now serviced by machines 
using electronic principles. In department stores, new controls supervise as 
many as O remote operations, tying in lighting systems, air-conditioning and 
a host of other units which previously required human guidance. 


The list of instances where automation is already used in one form or 
another is long. But the surface has hardly been scratched. In terms of the 
nation's total production, the output through automation is now but a 
pittance. 


Because automation and electronics are inseparable, it will definitely be 
a dominating factor in our own industry in the era ahead. Five years ago 
sales of automatic control devices by all companies amounted to a billion dol- 
lars..-In the past couple years, the figure has tripled, and it is still 
increasing. 


One recent survey showed more than 750 companies are now making coritrol 
systems, components and equipment. The entrance of so many new'firms into the 
field is predicated on the awareness that eventually all of industry must 
automate to some greater or lesser degree. There can be no doubt that many 
companies using self-regulating equipment will obtain a marketing edge over 
their competitors. The other must follow suit or fall by the wayside. When 
an employer does not automate,’ he faces bankruptcy and his workers face the 
loss of their jobs. 


In industry automatic production is based on the idea of unbroken motion 
and involves the combination of several piecemeal operations into one unit. 
By using the "building block" technique, you link a number of automatic sys- 
tems, gradually broadening the process until the completely automatic plant is 
within view. 


The ultimate, of course, would be the pushbutton type of operation cover- 
ing inspection, assembly, test, packaging and a raft of other services in 
which no human help is needed. We know of no completely automatic factory 
operating on a large-scale basis at the present time. The closest approach to 
it is the processing of some fluids -- in the chemical, oil and food indus- 
tries -- where only elaborate piping, pumping and time cycling is required. 
One leading organization of chemical producers has estimated that instruments 
and control mechanisms now commercially available can supply 90% of the needs 
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of the modern chemical plant. This is illustrative of the great things possi- 
ble with automation. But let us remember that while these things can be done, 
it doesn't mean they will be done. 


Complex machines to make improved goods quicker are costly. In the hur- 
ried pace of the future, marketing habits will change and some machines are 
bound to become obsolete. Management must look at the situation from the 
standpoint of the number of hours a specific machine can do the required job 
before it has to be altered. There was an instance of this in the automobile 
industry last fall. You may be familiar with the story of the car manufac-~ 
turer who went all-out for automation in producing engines -- to a certain 
point. Workmen still inserted semifinished connecting rods by hand into a 
device carrying them through transfer machines. Here's a case where availa- 
ble automatic equipment could have solved the particular problem. But at 
$90,000 it cost too much, and the manufacturer resorted to a combination of 
automatic and manual principles to get his finished product. 


On the other hand, while looking at the expense factor, there's the case 
of the company in the East which recently contracted for an electronic brain 
at a monthly cost of $30,000 -- $360,000 a year -- and figures it will be 
money well spent. In a matter of minutes, the machine supplies information 
that formerly required hundreds of manhours. The process used to take so long 
that figures were obsolete by the time they were computed. Now the system is 
used to speed the whole operation to the point where extra clerks are needed 
to keep business rolling once the desired answer is obtained. 


We contime to hear rumblings from people who feel that the machines of 
automation will bring contagious unemployment. They've compared it to a 
highly communicable virus that will have all of the nation's workers flat on 
their backs, when, in reality, the minor readjustments accompanying the early 
stages of automation should be considered as being sort of a vaccine. Karly 
immunization with the mild dose will bring lasting protection in the long run. 
In automation, we can look at the minor temporary dislocations as being but a 
step toward the stabilization of employment in the future. 


We have heard all kinds of suggestions from those who claim automation 
has its darker side. Using the theory that when machines replace men, the 
men who have no jobs cease to be customers and throw other men out of work, 
some people have even suggested that the government stop issuing patents on 
automatic devices. If such a step had been taken 15 or 20 years ago, when 
mechanization was growing by leaps and bounds, the workers now employed in the 
field of electronics might have been deprived of their present jobs. At GE 
alone, there are 5,000 people today who are making products that were unheard 
of 15 years ago. Their jobs exist only because we have employed the princi- 
ples of mechanization and automation as a means of lowering consumer costs. 


In face of this, we still have some people who insist automation is a 
sinister menace leading to economic catastrophe through disintegration of the 
labor market. Fortunately, however, the realization is now beginning to take 
hold in some quarters that the monster is nothing more than a Halloween 
spook -~ once you take off the false face, the person underneath isn't bad at 
all. To those who have fostered automation as a step to a better way of life, 
there is no surprise in recent pronouncements by some resistance leaders that 
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standards of living can be raised considerably if the problem is handled 
correctly. 


Take 1930 as a base year for the past quarter of a century of progress. 
Since then, production of machinery has increased about four times. Our gross 
national product has done the same. Retail prices have doubled, but at the 
Same time, average industrial earnings for employees have more than tripled. 
The number of hours worked each week has gone down and the number of people 
working in industry is now double the figure of 25 years ago. The portion of 
those doing hard, unskilled work has been cut in half. Better incomes have 
given the workers more purchasing power. Shorter hours have meant more pay 
for less work, and more leisure time. 


The whole concept of our new economy is that cheaper prices across the 
counter will enable the customer to satisfy basic needs, creating additional 
purchasing power. Past experience has shown that, as needs decrease, wants 
intensify. Knowing consumers will require more goods than we can now produce, 
our developmental laboratories must begin to anticipate the markets of the 
future, 


The potential is great. In the next seven years, the population of the 
United States may go as high as 18) million. The increase -~ 2); million -- 
gives us a tremendous market to cultivate. 


Despite the pessimistic outlook of some everything we have read in recent 
months points to a need for an augmented labor force in the years ahead, not a 
smaller one. Due to the low birth-rate of the 1930's the number of people 
entering the work force will be slim -- some of the experts have estimated it 
at only 11%. This, at a time when we must double production if we are to meet 
the requisites of the growing market. If we do not use automatic machinery, 
how else are we to achieve this goal? 


From this, you might conclude that we have nothing but clear sailing 
ahead. To convey this impression would be to present a picture that is not 
entirely true. You can produce as much as you want with automation, but even- 
tually the customer will get a yen for a different product. If he doesn't, 
it's up to you to provide the necessary stimulation. Otherwise, your market 
becomes saturated and you come to a dead end. 


The very fact that we must resort to new designs to maintain an even mar- 
ket should, in itself, serve as a handbrake against the quick emergence of the 
wholly automatic plant. In applying the advanced technology, the manufacturer 
must measure the pace of progress to determine how extensive the application 
should be. 


The key to widespread use of automation is flexibility. The ability to 
switch quickly from one type of assembly to another, with a minimum of down- 
time on equipment, is a prime requirement. This is especially true in the 
case of the small supplier, who is pressured into a race to keep up with the 
big producer's demand for new parts. We at GE have come a long way in this 
respect with our automatic component assembly system which is adaptable to the 
needs of the big producer and small manufacturer alike. 
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That's one of the faults we can find today. with almost every discussion 
taking place on the subject of automation generally. Too many people believe 
automation's broadest use. will be in the field where only the big boys play. 
But the "little fellow" is definitely in the picture. -In some cases, his 
problems will be more intense than those arising from the progressive process 
of the larger plants. In some instances, however, he will be spared many of 
the headaches. For example, he will not have to tie up large amounts of 
existing capital in research and creative engineering. That's a proposition 
the larger manufacturers will encounter and will have to measure up to as the 
use of automation widens. 


To achieve solid commercial success, those companies using new techniques 
in great degree must eventually look to the public for new money. If we are 
to entice the consumer to part with what he has saved through automation's s 
cost reductions, we must give them new designs to even out the purchase 
pattern. The sound conduct of business demands it. This is particularly true 
if fast-moving fields like electronics, where all the knowledge at our command 
must be sharpened to bring forth new designs to satisfy the appetite of the 
public for better living. 


This responsibility is. squarely on the shoulders of industry's engineers. 
Already, hundreds of ideas are coming from the drawing boards. There's one 
inherent danger, however. We must not only satisfy the whims of today's 
buyer -- the fellow whose tastes we know -- but of the potential purchaser of 
five to ten years from now. The public will dictate whether an item is obso- 
lete before it gets to the production stage. With our new economy, we could 
well be forced into a situation compelling more frequent introduction of new 
items. It means our engineers will have to pay closer attention to deadlines 
and performance dates in the future. Those who have been disrespectful of 
deadlines in the past may find there is no longer any time for head-scratch- 
ing. In getting the originality we seek quicker, there may be an extension of 
the present trend to "team engineering." 


With the strain for fresh ideas, we must assume that there will be a need 
for more engineers than we now have. For that reason, we may be rightfully 
disturbed by some of the statements made recently by those who have pinned the 
responsibility on management to provide the means of moving displaced produc- 
tion people to new jobs. While reaching for a greater share of the profits, 
they make no allowance whatsoever for the sizeable amounts of capital which 
must be diverted to engineering research to create the products which will 
make those new jobs available. Some industries now have a ratio of as many as 
100 engineers for every 1000 workers. There is every reason to believe that, 
with automation, this proportion will be changed considerably. 


While much of our emphasis has been on the laboratory engineer, we must 
be mindful, too, of the fact that the manufacturing engineer will take on 
somewhat of a different appearance as advancements are made. As automation 
comes in, there will be less of a need for specialization and a greater 
requirement for those trained in a number of engineering lines. We will 
require controls specialists to supervise entire assemblages of apparatus, so 
the complicated machines will remain in proper working order. They will have 
to know more than the intricate details of instruments and electronics and 
will be required, too, to know a good deal about general machine operations -- 
like welding, stamping, and other functions. 
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Much of this knowledge will be acquired on the job. But, in a greater 
measure , it must be supplied by our schools of technology. Systeme engineer- 
ing has opened a whole new field of technical education. 


I think the facts speak for themselves. We need not fear that which has 
so entranced the science fiction writers for some decades, "the invasion of 
the robots." We do have a need to recognize that the acceleration of techno- 
logy, which marches hand in hand with the concept of automation, places the 
engineer even more firmly on the management team. It places greater responsi- 
bilities upon his shoulders, for he must, as a member of that management team, 
help determine whether the product must ts designed for automatic production. 
The decision cannot be made on the basis of "Can it be done?" ‘but "When 
should it be done?" with full consideration for the human, the marketing, the 
financial and the profit problems involved. 


The responsibility also rests upon the engineer to help find a solution 
to the great and growing need for increased technical education and for find- 
ing the incentives that will attract qualified young people into technical 
fields. 


As a member of the management team, the engineer must help see to it that 
we do not become so wrapped up in the magic of the word "automation" that we 
ignore some of the areas that can provide equally important rewards to our 
economy and our standard of living. Management must find better methods. of 
handling marketing and distribution and human relations problems. 


89 


IRE ACTIVITIES IN THE FIELD OF AUTOMATIC CONTROL 


(As mentioned in the editorial, the IRE has contributed to the 
development of feedback control for a number of years, but only 
recently has it officially entered the automatic control field. 
This official activity, particularly in establishing control stan- 
dards, is described briefly by J. E. Ward, present chairman of IRE 
Subcommittee 26 on Feedback Control Systems. -- Editor's Note) 


IRE Technical Committee 26 on Feedback Control Systems was established in 1951 
under the chairmanship of Professor William M. Pease of the Massachusetts 
Institute of Technology. Among the initial objectives of the Committee were 
the establishment of standards on letter symbols and terminology for feedback 
control systems, methods of measurement and the promotion of an IRE profes- 
Sional group in this field. Most of these initial objectives have been 
attained. The IRE Professional Group on Automatic Control was formed in 195) 
under the chairmanship of Robert B. Wilcox, an original member of Technical 
Committce 26. The first two standards prepared by Committee 26 were approved 
in 1955 by the IRE Standards Committee. The first of these, "Graphical and 
Letter Symbols for Feedback Control Systems -- 1955," appeared in the Novem- 
ber, 1955, issue of Proceedings of the IRE. The second standard, "Terminology 
for Feedback Control Systems," appeared in the January, 1956, Proceedings. 


STANDARDS AND THE CONTROL FIELD 


In the IRE organization, the technical committees, working under the 
Standards Committee, are responsible for the preparation of all standards of 
terminology and methods of measurement. A committee cannot publish its work 
until completed and approved by the Standards Committee and the IRE Executive 
Committee, nor can it sponsor meetings. On. the other hand, a professional 
group cannot work on standardization, although it can sponsor meetings and 
suggest areas or material for consideration to the technical committees. The 
committees and the groups thus complement each other. 


AS one example of the need for standardization in the control field, con- 
Sider the name of the committee itself. Initially it was called the Technical 
Committee on Servo-Systems. Later, the name was changed to Feedback Control 
Systems, because this had become a more widely used term for the field. The 
parallel AIEE committee had made this same name change earlier. Whence comes 
the name "automatic control" used by the new professional group? This is an 
even broader term for our field and is also used as the name for ASME and 
American Standards Association committees in this field. This difficulty in 
even naming the field of interest indicates the need for adequate standards. 


The original committee consisted of six members, all of whom are still 
active. The current membership is seventeen. The Committee has at ‘the pres- 
ent time one Subcommittee 26.1 on Terminology for Feedback Control Systens, 
under the chairmanship of M. R. Aaron of the Bell Telephone Laboratories which 
has been working since 1952 on terminology. Mr. J. C. Lozier was the organ- 
izer and the first chairman of this subcommittee. A second subcommittee on 
methods of measurements is in the process of being formed. 
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The current committee membership is: 


J. E. Ward, Chairman _ Ee A. Sabin, Vice-Chairman 
M. R. Aaron DeeLerLippitt 

G. S. Axelby _ Je C. Lozier 

G. A. Biernson T. Kemp Maples 

R. &. Graham W. M. Pease 

V. B. Haas P. Travers 

R. Je Kochenburger R. B. Wilcox 

D. P. Lindorff So B. Williams 

W. Ke Linvill F. Re ZatLlin 


The current membership on Subcommittee 26.1 is: 


M. Re Aaron, Chairman 


G. S. Axelby T. Flynn 

cas ie Arthur dietitian Lozier 

V. Azgapetian C. F. Rehberg 
F. Zweig 


With the completion of the initial standards on symbols and terminology, 
the Committee is proceeding with the preparation of additional standards and 
material on standards of performance and methods of measurement. The Commit- 
tee feels that the expanding interest of the IRE in automatic control warrants 
intensive effort in these standardization activities. In addition to their 
activities on Committee 26, a number of committee members are also active in 
similar committees of other professional societies and in the American Stan- 
dards Association. This direct coordination of the work of committees of the 
various societies will hasten the realization of universal standards in the 
control field. 

John E. Ward 
Chairman, Committee 26.0 
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PROFESSIONAL GROUP ON AUTOMATIC CONTROL 


On October 5, 195, the IRE Board of Directors approved the petition for the 
formation of a Professional Group on Automatic Control. Sponsorship of this 
professional group was by members of the IRE Technical Committee on Feedback 
Control Systems and members of a number of IRE sections in answer to desires 
for a specialized organization which would be active in all phases of auto- 
matic control. A definite need was felt for an organization that would 
effectively aid in the tremendous progress being made in the field of auto- 
matic control by making contributions at the technical sessions of the IRE 
conventions, by organizing local professional chapter meetings, by publishing 
papers in the Transactions on Automatic Control, and by sponsoring national 
technical symposia. Ts . 


The purpose of the PGAC was given considerable thought in choosing the 
field of interest which appears in the constitution: 


"The field of interest of the Group shall be automatic control Sys- 
tems. It shall encompass the components thereof, such as transducers, 
data transmission links, computers and control devices and the integra- 

tion of these components into control systems. It shall include scien- 
tific, technical, industrial or other activities that contribute to this 
field, or utilize the techniques or products of this field, subject, as 
the art develops, to additions, subtractions or other modifications 
directed or approved by the IRE Committee on Professional Groups." 


The main aspect of the Groupts activity is the integration of computers, 
data processing, electronic techniques, feedback control, Switching circuits, 
telemetering, etc., all under the broad umbrella of automatic control. Such 
things as inventory control, manufacturing, accounting, and even business 
operations, can be considered control Situations in that they are amenable to 
automatic operations. Yet the scope has been maintained broad enough to cover 
the more general control problem with which the feedback control engineer is 
familiar. The real challenge is the integration of component parts and sub- 
systems into a complete control system. 


There are three main areas made available to the PGAC members (at present 
close to 1200) where they can obtain and disseminate information on the theory 
and practice of design of control systems: the national group symposia and 
technical sessions on a national or regional scale, the local professional 
chapters and the IRE Transactions on Automatic Control. 


fo date the PGAC has sponsored technical sessions at a number of conven- 
tions. At the 1955 National Convention the program included a technical | 
session on automatic control and the panel session entitled "Trends in Auto- 
matization of Procedures and Processes in Business and Industry," with Dr. G. 

- Brown of the Massachusetts Institute of Technology as chairman. The 

national group and the Dallas-Fort Worth Chapter sponsored the technical 
sessicn on automatic control at the South Western IRE Conference and Elec- 
tronies show in February, 1955. Again in August the national group with the 
Los Angeles Chapter sponsored the automatic control session at Wescon. These 
activities will become an annual event and will be supplemented by additional 
convention programs. When our Group becomes sufficiently strong, we plan to 
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initiate a symposium on automatic control, possibly in conjunction with one of 
our sister societies. 


Within the professional group organizational structure are the local 
professional chapters. The chapters have their own organization with a chair- 
man, vice-chairman, secretary-treasurer and activities committee, but maintain 
close coordination with the local sections. Such an organization has the def- 
inite advantage of providing local activity in one's professional field by 
means of meetings, panel sessions, etc., which are sponsored by the chapter or 
with other chapters in closely allied fields of interest. In this way, the 
members can take an active part in the group functions. A dozen or more 
active chapters will provide a strong national organization which in turn will 
be able to provide more for its members. During 1955, chapters have been 
organized in Baltimore, Boston, Dallas-Fort Worth, Los Angeles and Twin 
Cities. At present the New York Chapter is being organized as well as the 
Philadelphia, Connecticut Valley, San Diego and San Francisco. Monthly meet- 
ings have been arranged by these chapters to provide technical sessions for 
their members. Most of the papers presented to local organizations are also 
made available for publication in the Transactions on Automatic Control or the 
IRE Proceedings. 

\ 

The 1955 Convention Record Part IV contained those papers presented at 
the convention and represented the first published material of the PGAC. This 
spring 1956 issue of the Transactions on Automatic Control is our first Group 
publication. It is hoped that two or three such issues will be brought to the 
members during 1956 and on a quarterly basis thereafter. The success of the 
PGAC and its publication depends on the members to supply the Papers Procure- 
‘ment Committee with technical material. 


I would like to take this opportunity to inject a note of appreciation to 
the Administrative Committee, the Standing Committee members and Chapter 
Officers for their cooperation and effort since the formation of PGAC. A con- 
certed effort on the part of all members during the coming year will be neces- 
sary to increase our technical meeting activity and to provide quarterly 
issues of the Transactions. 

Robert B. Wilcox 
Chairman, PGAC 
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INFORMATION FOR AUTHORS 


George S. Axelby 
Westinghouse Electric Corporation 
Baltimore, Maryland 


Summary--Because the first PGAC Transactions papers will be reproduced from 
photographs of original manuscripts, 10 will be necessary to maintain a standard 
form to provide uniformity, neatness and legibility. 


This manuscript is written and typed in the desired form. Detailed instruc- 
tions for authors are included. 


INTRODUCTION 


Until the PGAC has grown sufficiently large, it will be necessary to publish 
the Group Transactions in the most inexpensive manner. This will be done by 
reproducing photographs of original, typewritten manuscripts and black-ink line 
“drawings. To achieve a desirable uniformity and legibility, the manuscripts 
should be typed with the same margins, spacing and headings wherever possible. 


Ordinarily a paper should contain the following sections: "Summary," "Intro- 
duction," main section (with appropriate headings), "Discussion," "Appendix" and 
"References." These will be considered individually in more detail. 


ARRANGEMENT OF MATERIAL 
Summary 


An abstract of no more than 200 words precedes the article and is called 
"summary." This section informs the reader briefly of the subject, method and 
results of the work. The summary should be informative and should outline the 
essentials of the paper. References or equations are not included. The essence 
of a good abstract is a concise statement of what was done and why. 


Introduction 


The purpose of the introduction is to orient the reader with respect to the 
problem. This section may be especially brief in a Transactions paper as much of 
the historical and introductory material will generally be known to Transactions 
subscribers, specialists in the field. 


Body of the Report 


In the main section or body of the report each step in the development of 
the material should be described as clearly and briefly as possible. Long math- 
ematical developments are to be placed in an appendix -- not in this section. 
Only a mathematical outline should be given and its importance and role in the 
development of the work clearly stated. 


Of course this section should not be headed "Body of the Report." But 
although it need not follow any particular form, for clarity there should be 
appropriate informative headings (and subheadings if necessary). 
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Discussion 


After the work has been described and the results have been given, the 
Significance of the work is discussed. In this part of the report some of the 
conclusions stated in the summary are developed. This section of the paper is 
important. Failure to state conclusions clearly is a common fault and the sig- 
nificance of a worthwhile contribution may be lost because the conclusions have 
not been clearly presented to the reader. 


Appendixes 
The appendixes include detailed subdevelopments of the main discussion. 


References 


A complete list of references should be given at the end of the paper and 
proper notation with superscripts made throughout the text. Each reference 
should be complete. For an article in a periodical the information contained 
should be: author's name, title of article, name of periodical, volume number, 
page numbers, month and year of publication. For a book: author's name, title of 
book, city of publication, publisher, year of publication, page numbers. 


Tllustrations and Tables 


Tables, photographs and diagrams are to be arranged so that they are self- 
explanatory and can be used with a minimum of reference to the text. The quanti- 
ties and units used in plotted curves should be given clearly and the legends 
worded to convey as much information as possible to the reader who may be casu- 
ally glancing through the Transactions. 


MARGINS, MEASUREMENTS AND STYLE 


All manuscripts are to be typed single space on anpelite typewriter (twelve 
characters to the inch) on 85 by 11 inch paper. The width of the text should be 
6 3/h inches exactly (81 characters) and the depth 8 3/4 inches. 


Titles and Subtitles 


The title of the paper should be concise, but clear. The author's name and 
affiliation are given after the title. The city and state of the author 's 
affiliation are also included, but street addresses and zone numbers omitted. No 
degrees or titles are mentioned. 


All main headings are also centered and capitalized, except for "summary." 
This particular heading is typed upper and lower case, underscored, and flush 


with the left-hand margin. A dash follows and the text begins on the same line. 
(See first page of this article.) 


Subheadings (or side headings) occupy a separate line, are underscored and 
upper and lower case. A line is skipped before and after the subheading. 


Tllustrations and Tables 


‘Al1 tables mst be typed and not drawn. Diagrams must be original drawings 
in black ink on white paper or tracing cloth; photographs mst be glossy prints. 
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Their size should not exceed 83 by 11 inches. Drawings with typewriting on them 
are not acceptable. All information mst be hand-lettered in black\ink, the 
lettering large enough to stand reduction. In graphs only the major coordinate 
lines should be drawn; a chart containing a large number of closely spaced lines 
will not reproduce well. 


Footnotes 


All footnotes are placed at the bottom of the page on which they are indi- 
cated. The first footnote is designated with an asterisk (+), the second with a 
dagger (+) and the third with a double dagger (+). 


SUBMITTING PAPERS 


The review and publication of papers will be seriously delayed if manu- 
scripts and illustrations are not submitted in proper form. The original and two 
copies should be sent to: 

Mr. M. R. Aaron 
Bell Telephone Laboratories 
Murray Hill, New Jersey 


Papers submitted should include a statement as to whether the material has been 
copyrighted, previously published or accepted for publication elsewhere. Mate- 
rial not acceptable for publication will be returned. 


DISCUSSION 
Using the information here and following the instructions as closely as pos- 


sible will result in neat, legible PGAC Transactions that should convey valuable 
information to the control engineer in a minimum amount of time. 
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